/W  fit 


(  '  '■ 


FINAL  REPORT 

(Grant  No.  AF0SR-82-0238) 


on 


DTIC 


ELECTE 
JUL  2  4  1986 

D 


HOT  ISOSTATIC  PRESSING  OF 
CERAMIC  POWDER  COMPACTS 


to 


AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 


February  7,  1986 


by 


J.  K.  McCoy  and  A.  J.  Markworth 


BATTELLE 

Columbus  Division 
505  King  Avenue 
Columbus,  Ohio  43201 


This  t&ohuioQi  r.v-vo  v  _ 

,or *• 
Distribution  UMira<Ji  IKVla. 

matthwj.  KMRna 

Chi.f,  Information  DiTi.nl on 


r_* vT\_* -rt*, n* rjA* -?\*~'v*^rxi<  *v* ajljul rx au^tx curj^^vf^jk a*.  -v* rjr^m ruwuv^> j\jp.b  jut* zuuti*jwww« 


February  13#  1986 


OBatreiie 

Columbus  Division 
SOS  Kint;  Avenue 
Columbus,  Ohio  020 1 
Telephone  ( (>  14  >  42l-(>424 
Telex  24-5454 


m 


*55 


p 

fe 


ha? 

HE 


*a 

v 

h 


S 


I 


•?; 


Major  Joseph  W.  Hager,  USAF 
Program  Manager 

Electronic  and  Materials  Sciences  Directorate 
AFOSR/NE 

Bolling  Air  Force  Base#  DC  20332 
Dear  Major  Hager: 


AF.Q.SR  Contract  SZ-Q223. 

Enclosed  are  sixteen  (16)  copies  of  the  final  report  on  the  above- 
mentioned  contract.  If  you  have  any  questions  regarding  the  report# 
please  feel  free  to  call  me  at  614-424-4050. 

Sincerely, 


J.  Kevin  McCoy 
Research  Scientist 
Physical  Metallurgy  Section 

JKM:taj 


Encs.  (16) 


N 


I 


j.^AVV*s - V SW JV .-VCVLy^V _N W. V .  <A!VxCV.V.xS%SVTVTg>2Vi  Cxr-4  • 


watiPi'  <  v.  **  'i  * ".'  fgw^mw  w.  yj  h»j  wiww  wklwj  w  vk 


TABLE  OF  CONTENTS 


Page 

I.  INTRODUCTION . 1 

II.  SUMMARY  OF  RESEARCH .  2 

III.  CONCLUSIONS .  6 

IV.  RECOMMENDED  RESEARCH  .  7 

V.  PRESENTATIONS .  8 

VI.  PUBLICATIONS  .  9 

APPENDIX  A 

A  RATIONAL  METHOD  FOR  CALCULATING  MECHANISM  MAPS  (reprint)  ...  A-l 

APPENDIX  B 

DILATOMETER  TECHNIQUE  FOR  HOT  ISOSTATIC  PRESSING  (reprir.)  ...  B-l 

APPENDIX  C 

CONTINUOUS  MONITORING  OF  VOLUMETRIC  CHANGES  IN  CERAMIC 

POWDER  COMPACTS  DURING  HOT  ISOSTATIC  PRESSING  (reprint) .  C-l 

APPENDIX  D 

INTERFACE-REACTION-CONTROLLED  KINETICS  IN  THE  HOT 

ISOSTATIC  PRESSING  OF  SUBMICROMETER  ALUMINA  POWDER  (reprint  .  .  D-l 

APPENDIX  E 

DISCUSSION  OF  "PRACTICAL  APPLICATIONS  OF  HOT-ISOSTATIC 

PRESSING  DIAGRAMS:  FOUR  CASE  STUDIES"  (reprint)  .  £-1 

APPENDIX  F 

DENSIFICATION  BY  INTERFACE-REACTION  CONTROLLED 

GRAIN-BOUNDARY  DIFFUSION  (preprint) .  F-l 

APPENDIX  G 

COMPUTER  SIMULATION  OF  EFFECTS  OF  THE  PORE  SIZE  DISTRIBUTION 

ON  THE  KINETICS  OF  PRESSURE-ASSISTED  FINAL-STAGE  DENSIFICATION 

OF  A  POROUS  SOLID  (preprint)  .  G-l 


security  cl 


PAGE  (Wh*t\  Dmtm  Entwd) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  REPORT  NUMBER 

AFOSR-TR-  86-0  348 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

3.  TYPE  OF  REPORT  9  PERIOD  COVEREO 

C Ol  OUN)  By- 

Final  C£c  XS 

4.  PERFORMING  ORG.  REPORT  NUMBER 

Hot  Isostatic  Pressing  of  Ceramic 

Powder  Compacts 

7.  AUTHORfA) 

«].  K.  McCoy  and  A.  J.  Markworth 

8.  CONTRACT  OR  GRANT  NUMBER^*; 

82-0238 

9. 


performing  organization  name  and  aooress 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


Battell e  Columbus  Laboratories 

505  King  Avenue,  Columbus,  Ohio  43201 

l^Os  (  ft 2- 

M.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

12.  REPORT  OATE 

January  15.  1986 

11.  NUMBER  OF  PAGES 

14.  MONITORING  AGENCY  NAME  a  ADDRESSf//  ditto  toot  from  Controlling  Oltlce) 

IS.  SECURITY  CLASS,  (ol  thle  report) 

UNCLASSIHED 

ISa.  DECL ASSIFICATION/ DOWNGRADING 
SCHEOULE 

16  DISTRIBUTION  STATEMENT  (ol  title  Report) 

W  ,  f  iX  "i\"  si 

Approved  for  publit  r»leaaa ; 
distribution  unlimited. 

IT.  DISTRIBUTION  STATEMENT  (ol  the  ebetrect  entered  In  Block  20,  II  dlllerent  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19.  KEY  WORDS  ( Contlnuo  on  rtvtric  old*  II  mcttiary  ind  Identity  by  block  number) 


Hot  Isostatic  Pressing,  Alumina,  solute  diffusion,  interface  controlled 


20.  ^ABSTRACT  (Continue  on  reverie  elde  If  '•erettory  end  .Identity  hy  hlnr.k  nurrt'T’> 

C  '  C 

he  densification  of  aluminum  oTraeTTn Trot  isostatic  pressing  has  been 
studied  in  detail.  Methods  for  emulating  maps  of  densification  rate 
as  function  of  temperature  and  applied  pressure  have  been  developed.  A 
new  mechanism,  interface-reacti^n-controlled  grain-boundary  diffusion, 
has  been  found  which  describes/the  densification  of  high-purity,  fine¬ 
grained  (grain  radius  of  0.7  pi)  aluminum  oxide  powder  at  temperatures 
up  to  1423  K.  Theoretical  models  have  been  developed  for  this  mechanism.. 


wnssro 


_ _  "V 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Entered}* 


DD  |  j2k*7J  1473  COITION  OF  I  NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAOEfWian  Data  E nf"4) 


for  both  the  initial  and  final  stages  of  densification.  Standard  geometries 
have  been  used  for  the  models;  the  initial  stage  is  described  in  terms  of 
impinging  spherical  particles,  and  the  final  stage  is  described  in  terms 
of  isolated  pores.  The  theoretical  model  is  found  to  describe  the  data 
for  fractional  densities  up  to  at  least  0.9,  although  it  breaks  down  for 
higher  densities.  From  theoretical  studies  of  final  stage,  we  find 
indications  that  the  discrepancy  between  theory  and  experiment  is  due  at 
least  in  part  to  neglect  of  the  effects  of  a  distribution  of  pore  sizes^ 


Accesion  For 

NTIS  CRA&I 
DTIC  TAB  . 
Unannounced 
Justification 


By  . 

Dist.  ibution  / 

Availability  Codes 

Dist 

Avail  ar-d/or 

Special 

M 

SECURITY  CLASSIFICATION  OF  THIS  PAGEfWian  Data  Enlaratf) 


K  AXj/Bt.lK&X'£s^i!iXXiX  "W.ryC  A^«fJtU\fLA%  PAZiAJ^nkAS  J1 


AAr^  k.  XJt  XIL^JSXan^iFTJk^  XT*. 


I 


1 


;/ 

ft 

i£s 


$ 

1  .% 

A1 


FINAL  REPORT 
on 

HOT  ISOSTATIC  PRESSING  OF  CERAMIC  POWDER  COMPACTS 
(Grant  No.  AF0SR-82-0238) 
to 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 

by 
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This  report  contains  a  summary  of  research  carried  out  under 
AFOSR  contract  No.  82-0238,  "Hot  lsostatic  Pressing  of  Ceramic  Powder 
Compacts",  the  period  of  performance  of  which  was  June  15,  1982  to  Janu¬ 
ary  15,  1986.  The  major  results  obtained  under  this  program  have  been, 
or  are  currently  being,  published  In  the  open  literature.  The  pertinent 
reprints  and  preprints  are  Included  as  appendices,  and  these  contain 
complete  descriptions  of  both  methods  and  results.  A  brief  summary  of 
the  research  is  presented  In  the  next  section,  with  references  made  to 
the  detailed  discussions  In  the  appendices. 
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SUMMARY .OF  RESEARCH 

The  objective  of  this  research  was  to  Investigate  the  poten¬ 
tial  of  theoretical  models  for  describing  and  predicting  the  denslflca- 
tlon  behavior  of  ceramics  In  hot  Isostatic  pressing  (HIP).  As  a  test 
material#  we  selected  Reynolds  RC-HP-DBM,  a  fine-grained,  high-purity 
alumina  powder.  The  kinetics  of  denslflcatlon  for  this  material  were 
studied  at  temperatures  of  1273  K  to  1423  K  and  pressures  of  34  MPa  to 
103  MPa.  The  primary  novel  feature  of  the  experiments  was  that  a  dlla- 
tometer  was  used  to  obtain  a  continuous  record  of  sample  density.  Use 
of  the  dllatometer  provided  a  very  detailed  view  of  the  kinetics  of  den¬ 
slflcatlon#  and  allowed  us  to  critically  evaluate  theoretical  models  for 
denslflcatlon. 

At  the  outset  of  this  work  we  expected  that  we  would  be  per¬ 
forming  a  large  number  of  simple  hot  Isostatic  pressing  experiments  and 
correlating  them  with  a  number  of  theories  of  denslflcatlon.  Accord¬ 
ingly,  we  developed  a  system  of  computer  programs  for  the  calculation  of 
denslflcatlon  mechanism  maps.  These  programs  are  self-contained  and 
Intended  to  be  readily  adapted  to  various  denslflcatlon  mechanisms  or 
even  to  different  problems  such  as  creep.  The  programs  are  also 
designed  to  be  "portable",  that  Is,  able  to  be  run  readily  on  a  variety 
of  computer  systems.  It  Is  expected  that  the  programs  would  run  on  many 
computer  systems  with  changes  to  only  the  few  routines  that  form  the 
Interface  with  the  plotter.  A  description  of  the  algorithms  used  In  the 
program  Is  found  In  Appendix  A. 

Study  of  the  denslflcatlon  maps  for  alumina  Indicated  that, 
for  applied  gas  pressures  of  up  to  100  MPa  and  temperatures  up  to  at 
least  1473  K,  the  dominant  mechanism  of  denslflcatlon  would  be  grain¬ 
boundary  diffusion  for  the  material  under  study.  Specifically,  the 
predictions  based  on  standard  theories  of  denslflcatlon  were  that  tran¬ 
sport  of  aluminum  Ions  would  be  rate-limiting,  that  aluminum  Ions  would 
move  primarily  by  grain-boundary  diffusion,  and  that  charge  compensation 
would  occur  primarily  by  diffusion  of  oxygen  Ions  through  the  lattice. 
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Near  the  end  of  the  first  year  of  the  project*  development  of 
the  HIP  dllatometer  was  completed.  This  offered  us  a  choice  of  two 
experimental  strategies:  we  could  perform  a  large  number  of  relatively 
inexpensive*  un Instrumented  HIP  experiments*  or  we  could  perform  a 
smaller  number  of  more  expensive  experiments  using  the  dllatometer  to 
monitor  changes  In  sample  geometry*  and  keeping  continuous  records  of 
temperature  and  applied  pressure.  A  few  tests  made  It  clear  that  the 
latter  choice  was  more  effective  due  to  the  large  volume  of  data  which 
could  be  collected  from  a  single  experiment.  In  addition,  when  taking 
data  with  the  dllatometer,  we  could  study  the  affects  of  changing  HIP 
conditions  without  considering  sample-to-sample  and  run-to-run  varia¬ 
tions.  A  discussion  of  the  principles  of  the  dllatometer  and  the  bene¬ 
fits  of  Its  use  Is  found  In  Appendix  B. 

Methods  of  data  analysis  were  developed  for  the  dllatometer. 
The  dllatometer  measures  the  diameter  of  an  encapsulated  cylindrical 
sample,  but  the  fractional  density  of  the  sample  Is  desired.  It  was 
necessary  to  develop  a  model  for  geometrical  changes  In  the  sample  and 
capsule  occurring  during  denslflcation  and  their  effect  on  the  measured 
diameter.  The  effects  of  denslflcation,  thermal  expansion,  and  changes 
In  the  length-to-dlameter  ratio  of  the  sample  were  Included.  Details  of 
the  data  analysis  are  described  In  Appendix  C. 

As  mentioned  above,  a  HIP  experiment  using  the  dllatometer 
provides  a  great  deal  of  Information,  while  an  ordinary  HIP  experiment 
gives  only  the  Initial  and  final  density.  In  particular,  the  dllatome¬ 
ter  provides  a  much  clearer  and  more  detailed  description  of  the  kinet¬ 
ics  of  HIP  than  Is  ordinarily  available.  To  accommodate  this  new  look 
at  HIP,  we  found  that  It  was  more  appropriate  to  consider  theoretically 
calculated  plots  of  density  as  a  function  of  time  than  It  was  to  use 
denslflcation  maps.  Accordingly,  we  developed  programs  to  solve  the 
differential  eauation  for  denslflcation.  This  approach  Is  complementary 
to  the  use  of  maps:  the  maps  provide  qualitative  Information  regarding 
changes  of  denslflcation  rate  In  response  to  temperature  and  pressure, 
while  Integration  of  the  denslflcation  rate  gives  the  details  of  the 
density  history  resulting  from  a  specific  HIP  treatment.  The  programs 
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for  calculating  density  histories  were  designed  so  that  experimental 
records  of  pressure  and  temperature  could  be  used  with  theoretical 
models  of  denslflcatlon.  As  a  result*  we  could  readily  plot  and  compare 
the  experimental  and  theoretical  denslflcatlon  histories.  This  provided 
a  very  graphic  means  for  evaluating  the  theory. 

The  most  Instructive  of  these  calculations  was  for  experiments 
In  which  we  changed  the  pressure  during  denslflcatlon.  In  standard 
theories  of  denslflcatlon  by  diffusion,  the  denslflcatlon  rate  varies 
linearly  with  the  applied  gas  pressure.  As  discussed  In  Appendix  D,  we 
found  that  our  data  could  not  be  reconciled  with  such  a  model.  However, 
the  data  were  consistent  with  a  model  In  which  the  denslflcatlon  rate 
varies  approximately  with  the  square  of  the  applied  gas  pressure.  The 
denslflcatlon  maps  calculated  early  in  the  project  were  Incomplete  In 
that  no  mechanism  with  such  a  pressure  dependence  had  been  Included. 

This  discovery  led  us  to  review  the  theoietlcal  basis  of 
models  for  a  wide  variety  of  denslflcatlon  mechanisms.  As  part  of  this 
review,  we  uncovered  an  error  In  a  published  model  of  denslflcatlon.  A 
description  of  the  correction  of  the  6rror  Is  found  in  Appendix  E. 

The  second-power  dependence  of  the  denslflcatlon  rate  on 
applied  pressure  forms  an  unusual  Intermediate  between  the  models  of 
denslflcatlon  by  diffusion,  which  predict  a  first-power  dependence,  and 
models  of  denslflcatlon  by  dislocation  motion,  which  predict  third-  or 
higher-power  dependence.  We  found,  however,  that  Interface-react Ion- 
controlled  grain-boundary  diffusion  had  been  proposed  as  a  mechanism  for 
creep  deformation,  and  that  this  mechanism  resulted  in  a  second-power 
dependence  of  creep  rate  on  applied  stress.  By  changing  the  grain 
geometry,  we  were  able  to  adapt  the  model  to  describe  denslflcatlon 
rather  than  creep.  Two  versions  of  the  model  were  developed:  one  with 
Impinging  spherical  particles  for  the  Initial  stage  of  denslflcatlon  and 
one  with  Isolated  spherical  pores  for  the  final  stage.  A  complete 
development  of  the  model  and  comparisons  of  e> ;on1mental  and  theoretical 
results  are  found  In  Appendix  F. 

The  model  of  denslflcatlon  by  interface-reaction-controlled 
grain-boundary  diffusion  was  found  to  fit  the  data  well  for  fractional 
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densities  up  to  O.S,  but  the  denslflcatlon  rate  was  overestimated  at 
higher  densities.  Since  the  final  stage  of  denslflcatlon  Is  important 
in  determining  the  mechanical  properties  of  high-performance  ceramics* 
this  situation  was  thought  to  deserve  closer  study.  We  found  that  both 
our  model  and  standard  models  of  denslflcatlon  by  diffusion  predict  that 
denslflcatlon  will  proceed  rapidly  to  completion.  This  Is  In  marked 
contrast  to  our  experimental  data  showing  a  gradual  decrease  In  denslfl¬ 
catlon  rate  and  to  the  common  observation  that  It  Is  often  difficult  to 
achieve  full  density. 

One  possible  reason  for  the  discrepancy  between  theory  and 
experiment  Is  that  standard  theories  assume  a  single  pore  size.  Since 
It  Is  assumed  that  all  the  pores  start  at  the  same  size  at  the  beginning 
the  final  stage  of  denslflcatlon*  they  all  remain  the  same  size,  shrink 
together,  and  finally  vanish  simultaneously.  This  Is  hardly  a  plausible 
picture.  A  more  reasonable  approach  Is  to  assume  that  there  Is  a  dis¬ 
tribution  of  pors  sizes,  with  small  pores  vanishing  rather  quickly  while 
large  pores  are  more  persistent.  Such  a  description  of  final-stage  den- 
slflcation  Is  given  In  Appendix  6.  It  Is  found  that  even  a  distribution 
of  pore  sizes  with  a  modest  variance  can  produce  a  much  more  gradual 
approach  toward  full  density. 
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The  following  conclusions  may  be  drawn  from  this  research: 

The  use  of  a  dllatometer  can  provide  a  detailed  record  of  den¬ 
sity  during  a  HIP  experiment.  This  amount  of  Information  obtained  In 
this  way  Is  so  much  larger  than  that  produced  from  an  ordinary  HIP 
experiment  that  It  Is  possible  to  reduce  greatly  the  number  of  experi¬ 
ments  and  total  cost  to  Identify  denslflcation  mechanisms  or  character¬ 
ize  denslflcatlon  behavior. 

While  denslflcatlon  maps  provide  a  good  summary  of  a  large 
body  of  denslflcatlon  data*  they  are  less  useful  for  analysis  of  experi¬ 
mental  kinetic  data.  With  data  from  a  dllatometer*  It  Is  more  appropri¬ 
ate  to  compare  plots  of  the  recorded  density  history  with  a  theoreti¬ 
cally  calculated  density  history. 

It  has  been  shown  that  models  using  a  single  particle  size  can 
provide  an  adequate  description  of  denslflcatlon  at  fractional  densities 
up  to  0.9.  At  higher  densities,  models  using  single  pore  and  particle 
sizes  predict  excessive  denslflcatlon  rates.  A  detailed  description  of 
the  microstructure,  especially  of  the  pore  size  distribution.  Is  neces¬ 
sary  for  adequate  predictions  of  final-stage  denslflcatlon  kinetics. 


RECOMMENDED  .RESEARCH 


As  mentioned  above,  models  of  denslficatlon  which  use  a  single 
grain  size  and  pore  size  do  not  provide  an  adequate  description  of  den- 
slflcatlon  at  fractional  densities  above  0.9.  Unfortunately,  the 
mechanical  properties  of  high-performance  ceramics  depend  sensitively  on 
density  In  this  range.  In  order  to  Improve  the  predictive  capabilities 
of  theoretical  models  to  the  point  at  which  they  can  be  used  effectively 
to  reduce  the  cost  of  developing  new  ceramic  materials.  It  is  recom¬ 
mended  that  additional  study  of  the  final  stage  of  denslficatlon  be  per¬ 
formed.  Since  It  has  been  found  that  the  final  stage  of  denslflcation 
is  si inif Icantly  affected  by  the  distributor  of  pore  sizes,  it  Is 
appropriate  to  perforin  detailed  quantitative  metallography  to  determine 
the  pore  size  distribution  functions  and  to  attempt  to  relate  these 
functions  to  powder  properties  such  as  particle  size  distribution  and 
degree  of  agglomeration. 
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PRESENTATIONS 

The  following  presentations  have  been  given  as  part  of  work  cn 
this  project: 

J.  K.  McCoy#  WA  Coherent  Method  for  Calculating  Mechanism  Maps”#  Basic 
Sciences  Division  of  the  American  Ceramic  Society#  Fall  Meeting# 
Columbus#  Ohio#  October  30-November  2#  1983. 

R.  R.  Wills  and  J.  K.  McCoy#  "Hot  Isostatic  Pressing  Denslflcatlon 
Maps"#  Basic  Sciences  Division  of  the  American  Ceramic  Society#  Fall 
Meeting#  Columbus#  Ohio#  October  30-November  2#  1983. 

R.  R.  Wi'ils  and  J.  K.  McCoy#  "Hot  Isostatic  Pressing  of  Al^  Powder: 
Theory  vs.  Experiment"#  Basic  Sciences  Division  of  the  American  Ceramic 
Society#  Spring  Meeting#  Pittsburgh#  Pennsylvania#  April  29-May  3#  1984. 

R.  R.  Wills,  L.  E.  Muttart#  and  J.  K.  McCoy#  "Continuous  Monitoring  of 
Denslflcatlon  and  Phase  Changes  During  Hot  Isostatic  Pressing"#  Basic 
Sciences  Division  of  the  American  Ceramic  Society#  Fall  Meeting#  San 
Francisco,  California#  October  28-31#  1984. 

J.  K.  McCoy  and  R.  R.  Wills#  "Dllatometer  Technique  for  Hot  Isostatic 
Pressing",  Ninth  Annual  Conference  on  Composites  and  Advanced  Materials# 
Cocoa  Beach#  Florida#  January  20-23#  1985. 

R.  R.  Wills  and  J.  K.  McCoy#  "Interface-Controlled  Kinetics  In  the  Hot 
Isostatic  Pressing  of  Alumina  Powder",  American  Ceramic  Society#  87th 
Annual  Meeting#  Cincinnati#  Ohio#  May  5-9#  1985. 
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INTRODUCTION 


Since  the  introduction  of  mechanism  maps  by  Weertman  and  Weertman  (!)  and  their  subse¬ 
quent  development  by  Ashby  (2),  there  has  been  a  great  and  growing  interest  in  the  construction 
and  use  of  such  maps.  They  have  been  used  to  describe  the  kinetics  of  creep  (2,3),  hot  pressing 
(4),  sintering  (5),  and  fracture  (6).  In  addition,  potential-pH  diagrams  (7),  which  are  con¬ 
structed  in  a  similar  way,  have  been  used  to  describe  the  corrosion  behavior  of  a  large  number 
of  materials.  The  discussion  in  tnis  paper  is  couched  in  terms  of  deformation  modeling,  but  the 
methods  may  be  applied  readily  to  other  problems. 

A  variety  of  algorithms  have  been  used  in  the  construction  of  these  maps,  but  for  the 
present  discussion,  it  will  suffice  to  classify  them  as  "linear"  or  "point-by-point".  The 
linear  methods  rely  on  the  use  of  straight  lines  to  represent  the  contours  and  boundaries  in 
the  map,  with,  in  some  cases,  the  use  of  a  change  of  variables  to  produce  the  necessary  linear¬ 
ization  (8).  The  point-by-point  methods,  on  the  other  hand,  avoid  the  assumption  of  linearity 
by  making  an  exhaustive  scan  of  the  whole  map.  The  point-by-point  approach  is  often  thought  of 
as  being  inefficient  because  past  programs  required  the  evaluation  of  rates  at  a  large  number  of 
points,  while  linear  methods  involve  the  treatment  of  a  linear  system  at  relatively  few  points. 
For  metallurgical  processes,  however,  it  is  noc  clear  that  linearization  is  always  possible,  so 
it  is  natural  to  choose  the  more  general  point-by-point  approach. 

Unfortunately,  the  point-by-point  programs  which  were  available  to  us  were  less  than 
inspiring.  They  did  not  meet  modern  standards  of  modularity,  documentation,  or  user  convenience, 
and  they  appeared  to  be  slipshod  in  their  use  of  data.  A  review  of  the  literature  further  sug¬ 
gested  that  no  study  had  been  made  of  the  methods  necessary  for  systematic  calculation  and  plot¬ 
ting  of  mechanism  maps.  Accordingly,  we  set  out  to  develop  methods  which  make  careful  use  of 
all  data  and  which  could  be  used  to  plot  naps  in  essentially  final  form.  Besides  the  intended 
advantage  of  producing  finished  maps,  we  also  found  that  our  algorithm  is  surprisingly  efficient 
and  requires  far  fewer  points  than  other  point-by-point  methods. 

Calculation  of  Boundaries 

Let  us  state  the  problem  of  the  construction  of  a  mechanism  map  in  abstract  terms.  We 
have  n  functions  f^,  i=l,  n,  each  of  which  is  dependent  on  two  variables.  At  any  point  (x,y), 
the  value  f-j(x,y)  is  the  rate  of  deformation  via  the  ith  mechanism.  The  map  describes  a  rec¬ 
tangular-  region  xm-n  <_  x  <_  x^x,  ymin  ^y  <_  ymx.  The  first  step  in  the  construction  of  the  map 
is  to  divide  the  map  into  one  or  more  regions  where  f-j(x,y)>  fj(x,y)  for  all  j  f  i.  Such  a 
region  will  be  called  a  region  of  dominance  of  mechanism  i.  Tne  second  step  is  to  find  a  set 
of  contour  lines  for  each  of  the  regions.  The  third  step  is  to  plot  the  boundaries  of  the 
regions  and  the  contours  in  a  smooth  and  efficient  manner.  In  the  remainder  of  this  paper  we 
will  discuss  each  of  these  tnrae  steps  in  sequence.  Finally,  we  will  present  examples  which 
illustrate  the  efficiency  of  the  algorithm. 

To  make  the  problem  tractable,  a  regular  grid  is  superimposed  on  the  map,  dividing 
the  area  of  the  map  into  rectangles.  The  density  of  the  grid  will  vary  from  one  map  to  another, 
based  on  the  complexity  of  the  map,  the  desired  level  of  smoothness  of  the  curves,  and  the  need 
for  computational  speed. 

In  calculating  the  map,  each  rectangle  is  treated  separately,  and  the  portion  of  any 
boundary  or  contour  which  lies  within  the  rectangle  is  represented  by  a  line  segment.  The 
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general  efficiency  of  the  algorithm  may  be  Improved  by  careful  storage  of  the  data.  It  is  de¬ 
sirable  to  evaluate  the  functions  fj  only  once  at  each  grid  point  and  to  save  the  results  of  the 
analysis  of  Individual  sides  of  a  rectangle  for  subsequent  use  in  analysis  of  adjacent  rectangles. 

Since  each  rectangle  is  treated  separately,  let  us  consider  the  determination  of  the 
boundaries  of  the  regions  of  dominance  which  fall  within  one  rectangle.  We  will  first  discuss 
the  topology  of  the  boundaries,  then  investigate  how  the  topology  and  positions  of  the  boundaries 
may  be  calculated. 

Figure  1  shows  four  topologically  distinct  arrangements  of  boundaries  between  regions 
of  dominance  in  one  rectangle  of  the  map.  We  have  found  that  real  mechanism  maps  very  rarely 
require  more  complex  topologies. 

In  calculating  the  positions  of  the  boundaries,  no  data  are  used  except  the  values  of 
the  functions  f^  at  the  four  corners  of  the  rectangle.  All  boundaries  within  the  rectangle  are 
assumed  to  be  line  segments. 

Rather  than  analyzing  the  rectangle  as  a  whole,  we  study  each  side  individually,  then 
combine  the  results  for  the  four  sides.  We  assume  that  the  rates  of  the  various  mechanisms  vary 
linearly  along  the  sides  of  the  rectangle.  Therefore,  if  one  mechanism  is  dominant  at  both  ends 
of  a  side,  it  is  dominant  over  the  entire  side.  If  two  different  mechanisms  are  dominant  at  the 
two  ends,  there  are  two  or  three  dominant  mechanisms  on  the  side,  and  the  points  where  the 
boundary  or  boundaries  Intersect  the  side  may  be  calculated  using  a  linear  interpolant.  Sides 
with  four  or  more  dominant  mechanisms  can  be  ignored  because  of  the  assumption  that  only  the 
four  topologies  of  Figure  1  are  present. 

The  topology  of  the  rectangle  may  be  determined  by  simply  adding  the  number  of  points 
where  boundaries  intersect  the  four  sides  of  the  rectangle.  Plotting  is  trivial  if  there  are 
zero  or  two  such  points.  If  the  boundaries  intersect  the  sides  of' the  rectangle  at  four  points 
(two  nonintersecting  boundaries),  it  is  only  necessary  to  draw  the  boundaries  between  appropriate 
pairs  of  points.  If  the  boundaries  intersect-the  sides  of  the  rectangle  at  three  points  (three 
intersecting  boundaries),  plotting  cannot  be  done  until  the  location  of  the  intersection  is 
found.  In  this  case,  we  temporarily  treat  the  rectangle  as  if  only  two  dominant  mechanisms  were 
present,  so  that  the  boundary  goes  all  the  way  across  the  rectangle.  The  three  line  segments  ob¬ 
tained  will  intersect  at  three  points,  the  centroid  of  which  is  used  as  the  intersection  of  the 
three  boundaries.  This  procedure  is  illustrated  in  Figure  2.  The  point  of  intersection  could  be 
calculated  directly  by  using  a  bilinear  interpolant  and  solving  for  the  point  at  which  all  three 
mechanisms  have  equal  rates. 

Calculation  of  Contours 

Let  us  now  consider  the  methods  used  in  calculating  the  positions  of  the  contours  of 
constant  deformation  rate.  First  we  will  consider  the  case  where  the  entire  rectangle  lies 
within  a  single  region  of  dominance,  then  we  will  see  how  the  procedure  must  be  modified  to 
treat  cases  in  which  the  rectangle  is  divided  between  two  or  more  regions  of  dominance. 

In  either  case,  the  first  step  in  contouring  is  to  find  out  which  contours  will  pass 
through  the  rectangle.  To  do  this,  we  must  find  the  minimum  and  maximum  rates  of  the  dominant 
mechanism(s).  If  a  single  mechanism  is  dominant  throughout  the  rectangle,  it  is  only  necessary 
to  check  for  extrema  at  the  four  corners.  If  there  are  two  or  more  dominant  mechanisms,  it  is 
also  necessary  to  check  the  rates  at  the  points  where  the  boundaries  intersect  the  sides  of  the 
rectangle  and.  If  there  is  such  a  point,  where  three  boundaries  intersect.  Given  the  minimum 
and  maximum  rate,  the  minimum  and  maximum  contour  levels  follow  immediately. 

For  the  case  in  which  a  single  mechanism  is  dominant  throughout  the  rectangle,  we  first 
search  all  four  sides  of  the  rectangle  for  points  where  the  rate  of  the  dominant  mechanism  is  the 
same  as  a  contour  level.  There  will  be  either  two  or  four  such  points.  If  there  are  two  points, 
the  contour  is  taken  to  be  the  line  segment  which  connects  the  points.  If  there  are  four  points, 
a  more  or  less  arbitrary  criterion  must  be  used  to  decide  which  pairs  of  points  will  be  connected 
by  contours.  t 

For  the  case  in  which  two  or  more  mechanisms  are  dominant  in  different  portions  of  the 
rectangle,  each  mechanism  is  contoured  individually  as  if  it  were  dominant  throughout  the  rec¬ 
tangle.  The  contour  segments  are  then  "cut  to  length"  by  truncating  them  at  the  point  where 
they  enter  a  different  region  of  dominance. 
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Sorting  and  Plotting 

If  the  methods  described  above  are  applied  to  every  rectangle  in  the  map,  the 
boundaries  and  contours  will  be  completely  calculated  -  ,  at  least  in  principle,  ready  to  be 

plotted.  To  avoid  inefficient  use  of  the  plotter  and  t..  assure  a  smooth  appearance,  however,  it 
is  necessary  to  sort  the  line  segments. 

It  is  clear  that  standard  algorithms  for  sorting  a  list  of  numbers  cannot  be  used. 

All  of  the  efficient  sorting  algorithms  (Shell  sort  (9),  guick  sort  (10),  heap  sort  (11))  make 
use  of  the  fact  that,  if  two  numbers  are  chosen  at  random  from  the  list,  it  is  possible  to  tell 
which  of  the  two  will  precede  the  other  when  the  sort  has  been  completed.  With  a  list  of  line 
segments,  however,  it  is  usually  not  possible  to  tell  if  two  randomly  chcsen  line  segments  are 
even  part  of  the  same  curve. 

The  simplest  approach  to  this  problem  would  be  to  compare  the  position  of  the  last 
known  point  on  the  curve  with  each  endpoint  of  each  line  segment  in  the  list.  Segments  which 
have  been  plotted  would  be  marked  in  some  way  to  prevent  subsequent  replotting.-  Such  an 
approach  would  work,  but  it  would  be  very  slow  if  the  list  were  long. 

The  efficiency  of  the  basic  algorithm  described  above  can  be  improved  in  two  ways: 
First,  we  can  restrict  our  definition  of  a  curve,  and  second,  we  can  exploit  the  existing  order 
in  the  list  of  line  segments. 

Since  we  are  dealing  with  a  list  of  line  segments,  the  "curves"  which  we  have  been 
discussing  are  in  fact  polygonal  paths.  Usually,  each  side  of  these  polygonal  paths  corresponds 
to  a  single  line  segment.  Let  us  now  define  a  "directed  curve"  as  a  polygonal  path  in  •..'hich 
each  side  points  to  the  right  or  straight  up.  Figure  3  shows  a  curve  and  its  dissection  into 
directed  curves. 

From  the  definition  of  a  directed  cu1  ve:  it  follows  that  each  of  the  line  segments  in 
a  directed  curve  has  a  definite  starting  poinc  and  a  definite  ending  point.  The  x  coordinate 
of  the  starting  point  is  less  than  the  x  coordinate  of  the  ending  point,  and  if  these  are  equal, 
the  y  coordinate  o.~  the  starting  point  is  less  than  the  y  coordinate  of  the  ending  point.  Also, 
each  line  segment  in  a  directed  curve  except  the  first  and  the  last  has  a  unique  predecessor  and 
a  unique  successor. 

The  use  of  directed  curves  has  two  practical  consequences  for  the  sorting  and  plotting 
of  a  list  of  line  segments.  First,  as  Figure  3  makes  clear,  complicated  curves  may  be  broken  up 
into  several  directed  curves  which  are  plotted  separately.  This  is  not  desirable,  tut  it  has 
little  impact  upon  the  preparation  of  mechanism  maps  since  it  rarely  occurs.  Second,  when 
searching  for  a  successor  to  the  current  line  segment,  the  position  of  the  ending  point  of  the 
current  line  segment  need  only  be  compared  to  the  starting  poin+s  of  the  potential  successors; 
the  ending  points  may  be  ignored.  This  cuts  the  length  of  the  search  in  half. 

Let  us  now  investigate  how  the  search  time  can  be  cut  further  by  exploiting  the  exist¬ 
ing  order  in  the  list  of  line  segments.  When  we  calculate  a  map,  we  scan  the  map  in  columns  of 
rectangles,  working  from  the  left  edge  of  the  map  to  the  right  edge.  Within  each  column,  we 
scan  from  the  bottom  of  the  column  to  the  top.  Within  each  rectangle,  we  sort  the  list  of  line 
segments  according  to  the  y  coordinate  of  the  starting  point  of  the  segment.  Since  the  rec¬ 
tangles  do  not  overlap,  the  ordering  of  the  y  coordinate  extends  throughout  the  column.  The 
list  of  line  segments  is,  therefore,  highly  ordered:  it  is  divided  into  columns;  the  segments 
in  one  column  are  ordered  from  bottom  to  top;  each  segment  points  to  the  right  or  straight  up. 

From  this,  it  is  clear  that  the  successor  to  a  given  line  segment  must  be  either  in 
the  current  column  (the  column  which  contains  the  current  l-'ne  segment)  or  in  the  column  just 
to  the  right  of  the  current  column.  Therefore,  it  is  only  necessary  to  search  two  columns  when 
looking  for  the  successor  rather  than  the  whole  list.  If  no  successor  is  found  in  these  two 
columns,  the  curve  ends. 

Further  gains  in  efficiency  can  be  realized  by  using  the  fact  that  the  segments  are 
ordered  on  the  y  coordinate  of  their  starting  points.  When  searching  the  current  column,  the 
search  may  start  at  the  current  line  segment  and  stop  as  soon  as  the  y  coordinate  of  o  potential 
successor  is  so  large  or  so  small  that  segments  which  are  further  up  or  duwn  the  column  cannot 
possibly  join  up  with  the  current  segment.  A  similar  strategy  may  be  used  in  searching  the 
column  just  to  the  right  of  the  current  column. 
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Examples 

Figure  4  shows  two  mechanism  maps  which  were  calculated  and  plotted  by  the  methods 
described  here.  The  same  data  and  constitutive  equations  were  used  for  both  maps,  but  grids  of 
two  different  densities  were  used.  It  has  been  claimed  (3,8)  that  a  typical  map  requires  the 
evaluation  of  deformation  rates  at  4000  to  6000  points,  and  maps  have  been  published  with  at 
least  6900  evaluation  points  (12).  By  contrast.  Figure  4a  used  evaluations  at  441  points,  and 
irregularities  in  the  curves  are  barely  perceptible.  Only  169  evaluations  were  used  in  the 

preparation  of  F’.jre  4b.  The  resulting  irregularities  in  the  map  are  visible  but  hardly  ob¬ 

jectionable.  Clearly,  careful  use  of  the  data  results  in  a  dramatic  improvement  in  efficiency. 

Acknowledgment 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under 
Contract  AF0SR-82-0238. 

References 

1.  J.  Weertman  and  J.  P..  Weertman,  in  R.  W.  Cahn  (ed.)  Physical  Metallurgy,  North-Holland, 
Amsterdam,  p.  793  (1965). 

2.  M.  F.  Ashby,  Acta  Met.,  20,  887  (1972). 

3.  T.  G.  Langdon  and  F.  A.  Mohamed,  Mater.  Sci.  Eng.,  32,  103  (197r). 

4.  P.  A.  Urick  and  M.  A.  Notis,  J.  Am.  Ceram.  Soc.,  56,  570  (1973). 

5.  F.  B.  Swinkels  and  M.  F.  Ashby,  Acta  Met.,  29,  259  (1981). 

6.  C.  Gandhi  and  M.  F.  Ashby,  Acta  Met.,  27,  1565  (1979). 

7.  P.  B.  Linkson,  B.  D.  Phillips  and  C.  D.  Rowles,  Corros.  Sci.,  19,  613  (1979). 

8.  T.  G.  Langdon  and  F.  A.  Mohamed,  J.  Mater.  Sci.,  13,  1282  (1978). 

9.  D.  L.  Shell,  Coran.  ACM,  2,  No.  7,  30  (1959). 

10.  C.A.R.  Hoare,  Comm.  ACM,  4,  321  (1962). 

11.  J.W.J.  Williams,  Comm,  ACM,  7,  347  (1964). 

12.  A.  Mohan,  N.  C.  Soni  and  V.  K.  Moorthy,  J.  Nucl .  Mater.,  79,  312  (1979). 


FIG.  1 

Four  Possible  Topologies  for  Boundaries  in  One 
Rectangle  of  the  Map 

a.  No  boundaries,  b.  One  boundary, 

c.  Two  nonintersecting  boundaries, 

d.  Three  intersecting  boundaries. 
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Treatment  of  Three  Intersecting  Boundaries:  a,b,c.  Location  of  single  boundary 
if  third  mechanism  were  absent.  Dotted  poition  of  boundarv  will  be  cut  off  by 
third  mechanism,  d.  Superposition  of  a,  b,  and  c.  Due  to  approximations, 
boundaries  do  not  intersect  at  a  single  point.;  e.  Composite  map  with  regions 
of  dominance  labeled.  Intersection  is  drawn  at  centroid  of  middle  triangle  in  d. 
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FIG.  3 

A  Typical  Curve  and  Its  Dissection  Into  Three 
Directed  Curves 
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Dilatometer  Technique  for  Hot-Isostatic  Pressing* 


J.  Kevin  McCoy  and  Roger  R.  Wills* 

Battelle  Columbus  Laboratories 

505  King  Avenue,  Columbus,  OH  43201-2693 


A  dilatometer  technique  is  described  for  hot  isostatic  pressing  (HIP)  that  provides 
far  more  data  from  each  HIP  experiment  than  can  be  obtained  by  conventional  HIP 
practice.  The  design,  data  analysis,  and  typical  results  obtained  with  the  dilatome - 
ter  are  described. 
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Before  hot  isostatic  pressing  (HIP)  can  be  applied  in  the  production  of  a 
ceramic  component,  it  is  necessary  to  determine  an  appropriate  HIP  cycle 
for  developing  the  desired  properties,  such  as  a  given  level  of  density. 
Those  who  are  involved  in  cycle  development  for  HIP  know  that  this  is 
frequently  a  long,  difficult,  trial-and-error  process.  This  paper  describes  a 
dilatometer  technique  for  HIP  that  provides  far  more  data  from  each  HIP 
experiment  than  can  be  obtained  by  conventional  HIP  practice.  Because 
more  information  is  collected  from  each  experiment,  the  amount  of  work 
necessary  for  cycle  development  is  greatly  reduced.  The  dilatometer 
technique  was  developed  at  Battelle-Columbus  with  support  from  the  Air 
Force  Office  of  Scientific  Research.  To  explain  how  this  method  works, 
we  will  first  describe  the  design  of  the  dilatometer.  then  the  data  analy¬ 
sis.  and  finally  give  an  example  of  typical  results  obtained  by  this  meth¬ 
od. 

Procedure 

The  design  of  the  dilatometer  is  shown  schematically  in  Fig.  1 .  A  cy¬ 
lindrical  sample,  shown  end-on  in  the  figure,  is  heated  in  a  cylindrical 
furnace.  The  sample  temperature  is  monitored  by  a  thermocouple.  From 
the  cold  zone,  well  outside  the  furnace,  two  tungsten  rods  run  up  to  the 
sample.  The  rods  are  attached  to  pivots,  which  are  shown  as  triangles, 
and  are  held  in  contact  with  the  sample  by  a  spring,  which  is  not  shown. 
The  bottom  ends  of  the  tungsten  rods  are  attached  to  a  linear  variable 
differential  transformer  (LVDT).  The  entire  apparatus  is  enclosed  in  an 
autoclave. 

The  thermocouple  and  LVDT  provide  a  continuous  record  of  the 
temperature  and  diameter  of  the  sample  as  functions  of  time.  These  data 
are  subsequently  used  to  calculate  a  record  of  density  as  a  function  of 
time.  To  understand  how  they  are  used,  it  is  necessary  to  take  a  closer 
look  at  the  sample  geometry.  Figure  2  shows  a  cutaway  view  of  the  sam¬ 
ple.  The  sample  proper  is  a  ceramic  cylinder:  it  is  encased  in  a  close- 
fitting.  airtight,  cylindrical  metal  capsule.  Strictly  speaking,  the  dilatome¬ 
ter  provides  a  record  of  the  capsule  outside  diameter  rather  than  a 
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record  of  the  diameter  of  the  sample  itself.  Therefore,  the  computational 
procedure  may  be  summarized  as  follows:  We  assume  some  fractional 
density,  then  calculate  the  corresponding  capsule  diameter,  using  the 
treatment  described  below.  If  the  calculated  diameter  agrees  with  the 
value  reported  by  the  dilatometer.  then  the  assumed  fractional  density  is 
correct.  If  not.  a  new  fractional  density  is  chosen  and  the  process  is  re¬ 
peated. 

There  is  an  apparent  weakness  in  this  computational  procedure:  The 
dilatometer  gives  only  the  diameter  of  the  capsule,  but  the  length  will 
also  affect  the  calculated  density.  This  difficulty  is  handled  in  the  fol¬ 
lowing  way.  As  the  sample  densifies,  its  length  changes  from  its  initial 
value  L„  through  some  intermediate  value  L,  to  a  final  value  Lf. 

Similarly,  the  diameter  changes  from  its  initial  value  £),,  through 
some  intermediate  value  D,  to  a  final  value  Df.  The  initial  dimensions  can 
be  determined  before  the  sample  is  placed  in  the  capsule,  and  the  final 
dimensions  can  be  determined  after  HIP  by  removing  the  capsule.  It  is 
then  assumed  that  the  current  length  of  the  sample  is  related  to  the  cur¬ 
rent  diameter  by  the  equation 


The  value  of  the  constant  k  is  obtained  by  inserting  the  final  values  of  L 
and  D,  and  solving  for  k.  Given  Eq.  1  and  the  initial  density  and  initial  di¬ 
mensions  of  the  sample,  the  length  and  diameter  of  the  sample  can  be 
calculated  for  any  specified  density.  Typical  values  of  k  are  about  0.85.  so 
changes  in  the  length/diameter  ratio  are  relatively  small. 

Changes  in  the  geometry  of  the  sample  naturally  give  rise  to  changes 
in  the  geometry  of  the  capsule:  as  densification  proceeds,  there  is  a  dis¬ 
tinct  thickening  of  the  capsule  sidewall.  This  effect  is  treated  by  applying 
the  law  of  conservation  of  matter.  It  is  assumed  that  any  matter  in  the 
sidewall  at  the  beginning  of  HIP  stays  in  the  sidewall,  and  that  matter 
from  the  capsule  ends  does  not  transfer  to,  or  from,  the  sidewall.  It  is  also 
assumed  that  the  capsule  stays  in  close  contact  with  the  sample,  so  the 
inside  diameter  of  the  capsule  equals  the  diameter  of  the  sample,  and  the 
length  of  the  capsule  sidewall  equals  the  length  of  the  sample. 

The  computational  procedure  given  above  can  now  be  described  in 
greater  detail.  A  fractional  density  is  assumed  for  the  sample,  and  from 
this  the  sample  dimensions  at  20°C  are  calculated.  The  coefficient  of 
thermal  expansion  is  then  used  to  determine  the  sample  dimensions  at 
the  actual  HIP  temperature.  From  the  initial  sidewall  mass,  the  sidewall 
volume  at  20°  C  is  calculated,  and  the  coefficient  of  thermal  expansion  for 
the  capsule  material  is  used  to  determine  the  sidewall  volume  at  tempera¬ 
ture.  The  sample  dimensions  and  capsule  sidewall  volume  are  then  used 
to  calculate  the  capsule  sidewall  thickness,  and  the  capsule  outside  di¬ 
ameter  is  just  the  sample  diameter  plus  twice  the  sidewall  thickness.  If 
the  calculated  capsule  diameter  is  equal  to  the  diameter  reported  by  the 
dilatometer.  then  the  assumed  density  is  correct;  otherwise,  it  is  necessa¬ 
ry  to  try  a  different  fractional  density. 

Results 

The  type  of  results  obtained  by  this  method  is  shown  in  Fig.  3.  Dur¬ 
ing  most  of  the  heatup  of  the  sample,  essentially  no  densification  occurs. 
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While  the  sample  was  held  at  1050°C  and  34  MPa.  it  densified  to  about 
80%  of  its  theoretical  density.  At  this  point,  the  pressure  was  then  in¬ 
creased  to  100  MPa  in  order  to  speed  up  densification.  These  conditions 
were  maintained  until  the  sample  density  reached  91%  of  theoretical 
density.  Another  increase  in  densification  rate  was  then  obtained  by 
increasing  the  temperature  to  1 150°C.  Under  these  conditions,  the  sample 
reached  a  density  of  97%.  All  this  information  was  obtained  from  a  single 
experiment.  If  we  had  used  conventional  HIP  practice,  one  experiment 
would  have  given  us  only  the  initial  and  final  densities.  To  obtain  a  densi¬ 
ty  vs  time  curve  of  the  type  shown  here  would  have  required  many  exper¬ 
iments.  To  make  matters  worse,  there  would  have  been  sample-to-sample 
and  run-to-run  variations,  which  would  have  greatly  complicated  data 
analysis.  With  the  dilatometer.  however,  a  single  plot  is  obtained  which 
clearly  shows  the  response  of  the  system  to  changes  in  temperature  and 
pressure,  and  that  is  vital  information  in  developing  an  appropriate  HIP 
cycle. 

Conclusions 

In  summary,  the  primary  feature  of  this  dilatometer  technique  is  that 
it  provides  a  continuous  record  of  density  as  a  function  of  time  during 
HIP.  The  large  amount  of  information  provided  by  the  method  is  the  key 
to  its  principal  benefit:  much  faster  development  of  HIP  cycles. 

•A  more  detailed  version  of  this  work  has  been  accepted  for  publication  in  the  American 
Ceramic  Society  Bulletin. 

'Now  at  TRW,  Inc..  Cleveland.  OH. 
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The  hot  isostatic  pressing  (HIP)  of  ceramic  powder  compacts 
■  consists  of  several  sequential  steps.  These  are:  1)  green  preform 
preparation  by  an  appropriate  forming  process;  2)  preform  encap¬ 
sulation  tn  a  pressure-transmitting  membrane,  capsule,  or  “can"; 
3)  actual  HIPing  by  exposing  the  encapsulated  powder  compact  to 
a  high-pressure  gas  (typically  100  to  200  MPa)  at  elevated  tem¬ 
perature;  and  4)  removing  the  encapsulant  surrounding  the  HIPed 
powder  compact. 

Research  in  this  area  is  usually  performed  on  a  trial-and-error 
basis  by  analyzing  information  obtained  on  samples  before  and 
after  HIP'ing  in  an  effort  to  determine  the  optimum  HIP  conditions 
needed  to  achieve  either  complete  densification  or  to  produce  a 
ceramic  with  a  specific  fractional  density  (the  ratio  of  the  volume 
of  solid  in  a  material  to  the  total  volume  of  material  including 
pores).  The  process  is  iterative  and  sometimes  lengthy. 

The  technique  described  herein  allows  volumetric  changes  to  be 
monitored  continuously.  Data  are  collected  throughout  the  HIP 
cycle,  consequently  the  temperature  and  pressure  processing  pa¬ 
rameters  car  be  changed  if  the  selected  parameters  do  not  appear 
to  be  achieving  the  desired  objective  (usually  full  density).  Data 
output  is  comected  for  thermal  expansion  effects  to  obtain  exact 
density  changes  as  a  function  of  time.  This  paper  briefly  describes 
the  apparatus  and  the  error-correction  routine.  The  use  of  the 
equipment  is  demonstrated  in  HIP’ing  submicron  alumina  powder 
in  a  stainless  steel  capsule. 

Apparatus 

HIP  equipment  consists  of  gas  storage  tanks,  compressors,  and 
a  water-cooled  autoclave  containing  the  furnace  and  heat  shield 
package.  The  furnace  is  usually  loaded  through  the  top  and  the 
electrical  connections  lead  from  the  bottom  of  the  vessel.  A 
10. 16-cm  diameter.  1 1  68-cm  high  molybdenum  furnace  capable 
of  operating  at  1750°C  was  modified  to  contain  a  pedestal  and  a 
two-probe  dilatometer  (sec  Fig.  i)  attached  to  a  linear  variable 
differential  transformer  (LVDT).  The  3. 175-mm  diamea  i  tungsten 
probes  are  in  close  contact  with  a  stainless  steel  capsule  containing 


A  two-probe  diiatometer,  fitted  into  a  molyb¬ 
denum  furnace,  permits  continuous  record- 
inn  of  the  volumetric  changes  occurring 
inside  a  metal-encapsulated  ceramic  pow¬ 
der  compact  during  hot  isostatic  pressing. 
The  effects  of  thermal  expansion  of  the 
metal  capsule  and  the  ceramic  compact  are 
deleted  from  the  raw  data  by  a  computer 
program,  and  the  nonlinear  shrinkage  rela¬ 
tion  between  the  capsule  diameter  and 
length  is  allowed  for.  Use  of  the  equipment 
is  demonstrated  in  the  hot  isostatic  press¬ 
ing  of  a  submicrometer  alumina  powder  en¬ 
closed  in  a  stainless  steel  capsule. 


a  close-fitting  cylindrical  alumina  powder  compact.  Two  thermo¬ 
couples  monitor  the  temperature,  one  in  contact  with  the  platform 
supporting  the  specimen,  the  other  slightly  above  the  specimen. 

An  asbestos-based,  ceramic  insulating  plate  below  the  base 
heater  and  radiation  baffles  attached  to  the  tungsten  probes  main¬ 
tain  the  operating  temperature  of  the  LVDT  at  180°C  during 
HIPing  of  aiumina  at  1200°C.  The  LVDT  range  of  0.254  cm  is 
linear  to  within  0.25%.  Preliminary  experiments  were  performed 
using  a  molybdenum  disk  to  determine  the  accuracy  of  the  probe/ 
LVDT  equipment  under  HIP  conditions  A  multichannel  chart 
recorder  monitored  the  LVDT  output,  temperature,  and  pressure. 
The  measured  thermal  expansion  of  this  disk  at  1300°C  and 
1600°C  differed  from  the  published  data1  by  0.5%. 

Experimental 

Cylindrical  specimens,  26-mm  diameter  by  102  mm  long,  were 
formed  from  submicron  alumina*  powders  by  isostatic  pressing  at 
345  MPa.  These  specimens  were  then  bisque  fired  at  1120°C  and 
ground  to  final  22-mm  diameter  and  39.68  mm  length.  This  proce¬ 
dure  was  adopted  to  ensure  that  the  alumina  specimens  closely  fit 
the  stainless  steel  capsules.  The  specimens  and  caprules  were  then 
dried  and  vacuum  outgassed  at  1000°C  for  1  hour.  Each  specimen 
was  weighed,  and  the  initial  dimensions  of  each  specimen  and 
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Table  I.  Physical  Measurements  ot  Samples  and  Capsules  for 
Hot  Isostatic  Pressing  Required  for  Data  Manipulation 


Experiment  No. 

Dimension  (mm) 

i 

2 

3 

4 

5 

6 

Initial  sample*  length 

39.67 

39.68 

39.70 

39.66 

39.67 

39.68 

Initial  capsule  length 

51.32 

52.24 

49.90 

50.91 

50.82 

52.45 

Final  capsule  length 

46.82 

47.90 

45.28 

49.04 

46.35 

48.22 

Initial  total  end  plug 
length 

11.10 

11.13 

9.45 

10.97 

11.10 

12.76 

Final  total  end  plug 
length 

11.69 

11.60 

10.08 

11.16 

11.79 

13  57 

Initial  sample  diameter 

22.00 

22.00 

22.00 

22  00 

22.01 

22  01 

Initial  capsule  diameter 

25.39 

25.37 

25.38 

25.36 

25  37 

25.38 

Final  capsule  diameter 

23.28 

23.85 

23.33 

24.53 

23.11 

23.25 

Initial  capsule  wall 
thickness 

1.57 

1.57 

1.57 

1.57 

1.66 

1.66 

Final  capsule  wall 
thickness 

2.12 

1.98 

2.08 

1.83 

2.18 

2.72 

Initial  fractional 
density 

0.6562 

0.6556 

0.6556 

0.6549 

0  6370 

0.6368 

•Sample  refers  to  ceramic  compact 
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Fig.  1.  Schematic  of  diiatometer  and  furnace  for  hot  isostatic 
pressing.  Key:  1 .  encased  thermocouples  (W/Re6);  2.  specimen; 
3.  specimen  trough;  4.  thermal  barrier  and  cover;  5.  radiation  and 
gas  flow  baffles;  6.  perforated  moly  hearth;  7.  moly-tubular  pedestal; 
8.  moly  ribbon  heater;  9.  alumina  base  plates;  10.  radiation  shields; 
11.  couple  connection;  12.  power  connection;  13.  furnace  base; 
14.  transite  base;  15.  frame;  16.  tungsten  probes;  17.  inverted 
pedestal;  18.  spring  c*eel  hinge  each  end;  19.  extention  of  probe 
arms;  20.  very  light  compulsive  spnng;  21.  yoke;  22  spring  steel 
yoke  arms;  23.  physical  zero  adjust  screw;  24.  LVDT;  25.  core. 


capsule  taken  with  a  micrometer  Figure  2  shows  the  measured 
dimensions  and  actual  values  are  listed  in  Table  I  The  capsules 
were  then  sealed  under  vacuum  by  electron  beam  welding. 

After  each  sample  was  placed  in  contact  with  the  tungsten 
probes,  adjustments  were  made  to  the  LVDT  and  the  chart  recorder 
was  adjusted  to  zero.  The  zero  position  on  the  chart  recorder  was 
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Fig.  2.  Cross  section  of  capsule  and  sample  for  hot  isostatic  press¬ 
ing,  showing  measurements.  (1)  Sample  length,  (2)  capsule  length, 
(3,4)  end  plug  length,  (5)  sample  diameter,  (6)  capsule  diameter,  (7) 
capsule  wall  thickness.  Gap  between  capsule  and  sample  (shaded 
area)  is  exaggerated;  dashed  line  is  axis  of  capsule. 


set  to  accommodate  initial  contraction  of  the  capsule  diameter. 
Each  HIP  cycle  involved:  1)  initial  gas  pressurization,  2)  heating 
at  400CC  per  hour  to  the  desired  temperature  under  a  constant 
pressure,  3)  HIPing  under  the  required  temperature/pressure  con¬ 
ditions,  and  4)  allowing  the  specimen  to  cool  to  room  temperature. 
The  final  dimensions  of  the  capsule  (diameter  and  length)  and  the 
thickness  of  the  capsule  end  plugs  were  determined  with  a  micro¬ 
meter  (see  Table  I).  While  seven  measurements  are  necessary  to 
describe  the  initial  capsule  and  sample  geometries,  only  five  are 
required  after  HIPing  since  there  is  no  longer  a  gap  between  the 
sample  and  the  capsule.  The  immersion  technique  was  used  to 
determine  the  sample  density  after  the  capsule  was  removed. 

During  the  HIP  cycle,  changes  in  the  position  of  the  probe1:  are 
registered  in  the  LVDT  output  recorded  by  the  chart  recorder. 
Several  factors  are  responsible  for  this  movement.  Initial  pressur¬ 
ization  of  the  autoclaves  caused  some  small  movement,  but,  after 
pressure  equilibration,  the  chart  reading  returned  to  the  zero  posi¬ 
tion.  Dunng  heatup,  capsule  deformation  results  in  a  thickened 
capsule  wall  and  both  the  stainless  steel  capsule  and  the  alumina 
powder  compact  expand.  All  of  these  effects  give  rise  to  a  decrease 
in  the  value  registered  by  the  chart  recorder.  Densification  result* 
in  a  decreased  capsule  diameter,  and  the  chart  pen  moves  in  the 
opposite  direction.  The  actual  point  recorded  by  the  chart  pen  at 
any  one  time  is  the  sum  of  these  effects. 

Conversion  of  LVDT  Output  to  Fractional  Density  Data 
Data  from  several  sources  are  combined  to  produce  density  as  a 
function  of  time.  The  thermocouples  and  LVDT  provide  the  tem¬ 
perature  of  the  sample  and  the  outside  diameter  of  the  capsule  as 
a  function  of  time.  The  dimensions  of  the  sample  and  capsult  at  the 
beginning  and  end  of  the  HIP  cycle  are  also  needed  to  calculate  the 
(kr.fitv.  These  data  are  then  processed  by  a  program  that  compen¬ 
sates  for  thermal  expansion  and  thickening  of  the  capsule  wall. 

Data  reduction  *:?ns  by  keying  temperature,  pressure,  and 
capsule-shrinkage  data  ffon;  *he  strip  chart  recorders  into  a  com¬ 
puter.  Data  points  are  chosen  so  tnat  discrete  points  will  provide  a 
good  description  of  the  actual  experiment.  Since  the  time  interval 
between  data  points  is  small,  linear  interpolation  ;s  used  between 
the  readings.  Physical  dimensions  and  initial  density  dc*a  are  in¬ 
serted  at  the  beginning  of  the  file. 

The  data  are  converted  to  fractional  densities  by  a  FORTRAN 
com)  a  ter  program.  Before  discussing  the  algorithm,  it  is  appropri¬ 
ate  to  describe  the  underlying  assumptions  used  in  deriving  the 
algorithm  and  to  estimate  the  error  incurred.  Using  the  example  of 
an  alumina  powder  compact  encapsulated  in  a  stainless  steel  cap¬ 
sule,  the  integrated  thermal  expansions  of  alumina  and  stainless 
steel  are: 

aAlj03=7.2978£  -6(T-293) 

+7.8486E  — 1(XT-293)(T— 800)  (1A) 

c304SS=  1 ,7751£-5(T--293) 

+  5.0423£-9(T— 293KT-800)  (IB) 
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Table  II.  Capsule  Sidewall  Masses  Determined 
by  Various  Methods 


Experiment 

No 

Sidewall  mass  (g) 

Method  I 

Method  2 

Method  3 

1 

37.440 

36.929 

39.186 

2 

38.233 

36.903 

39.034 

3 

37.634 

36.932 

38.619 

4 

37.125 

36.872 

39.139 

5 

38.883 

38.843 

39.169 

6 

38.714 

38.704 

38.582 

In  these  equations,  which  are  integrals  of  the  ordinary  coeffi¬ 
cients  of  linear  thermal  expansion,  E  denotes  multiplication  by  the 
specified  power  of  10,  following  the  standard  computational  con¬ 
vention.  The  temperature  T  is  given  in  degrees  Kelvin.  These 
equations  were  obtained  by  fitting  a  quadratic  to  data  compiled  by 
Touloukian.'  Several  other  thermal  effects  might  be  anticipated. 
First,  the  probe-rod  dimensions  will  change  both  above  and  below 
the  pivots.  Second,  the  pivots  themselves  may  move  due  to  thermal 
expansion.  Third,  the  sample  may  be  displaced  as  its  supporting 
structure  expands.  Finally,  when  the  temperature  gradient  has  a 
component  normal  to  the  axis  of  the  probe  rods,  the  rods  will 
deflect.  It  is  not  feasible  to  perform  an  accurate  calculation  of  all 
thermal  expansion  effects.  Our  calculations  have  used  a  simple 
approach  by  assuming  that  all  thermal  expansions  except  those  of 
the  sample  and  the  capsule  have  a  negligible  net  effect  on  the 
measured  capsule  diameter.  The  validity  of  this  assumption  is 
supported  by  the  good  agreement  between  the  measured  thermal 
expansion  of  the  molybdenum  disk  and  the  literature  data. 

During  the  HIP  cycle,  the  thickness  of  the  capsule  sidewall 
changes  due  to  thermal  expansion  of  the  encapsulant  and  densi- 
fication  of  the  sample.  In  determining  sample  density,  it  is  neces¬ 
sary  to  consider  both  effects  and  to  make  some  assumptions  regard¬ 
ing  redistribution  of  matter  in  the  sidewall.  We  have  assumed  that 
the  sidewall  does  not  slip  axially  along  the  sample.  Thus,  changes 
m  the  length  or  diameter  of  the  sample  or  the  capsule  are  assumed 
to  be  taken  up  by  wall  thickening,  rather  than  by  moving  matter  to 
or  from  the  ends  of  the  capsule.  This  assumption  should  be  realistic 
since  the  capsule  is  always  being  forced  against  the  sample  by  the 
applied  gas  pressure.  Although  the  assumption  may  be  somewhat 
inaccurate  near  the  ends  of  die  sample,  our  measurements  of  cap¬ 
sule  diameter  are  made  near  the  middle  of  the  sample,  well  re¬ 
moved  from  end  effects. 

The  importance  of  the  correction  for  thermal  expansion  can  be 
seen  from  the  following  example.  Using  the  final  dimensions  from 
experiment  No.  1 ,  the  combined  thermal  expansion  of  the  sample 
and  the  capsule  upon  heating  from  20°  to  1150°C  is  0.378  mm  in 
diameter  if  the  fractional  density  of  the  sample  remains  constant. 
To  produce  this  large  a  change  in  size  at  constant  temperature 
would  require  a  change  of  0.08  in  the  fractional  density  of  the 
sample. 

Another  important  consideration  is  the  effect  of  temperature  and 
densification  on  the  shape  of  the  sample.  We  have  assumed  that  the 
thermal  expansion  of  alumina  is  isotropic  for  a  fixed  fractional 
density.  However,  the  length-to-diameter  ratio  of  the  sample  at  the 
beginning  of  the  HIP  cycle  is  different  from  that  at  the  end,  so  a 


Table  III.  Comparison  of  Final  Densities  as  Determined  by 
Calculation  from  Probe  Data  and  by  Direct  Measurement 


Experiment 

No 

Final  ft  actional  density 

Calculated 

Measured* 

Difference 

i 

0.9576 

0  960 

-0.002 

2 

0.8712 

0.876 

-0.005 

3 

0.9678 

0.952 

0.016 

4 

0.7697 

0.758 

0.012 

5 

0.9678 

0.984 

-0.016 

6 

0.9617 

0.970 

-0.008 

•Determined  by  liquid  immersion  technique 


change  in  shape  accompanies  the  density  change.  The  capsule's 
diameter  is  recorded  continuously,  but  its  length  is  measured  only 
at  the  beginning  and  end  of  the  HIP  cycle.  It  is  therefore  necessary 
to  assume  some  relationship  between  the  length  and  diameter  dur¬ 
ing  the  HIP  cycle.  It  is  assumed  that: 

L(p,20)  /  d(p,20)  \* 

I(p„20)  \d(p„20)J  { 

where  d(p, 20)  and  Z.(p,  20)  denote  the  sample  diameter  and  length, 
respectively,  at  the  fractional  density  in  question  and  corrected  to 
20°C.  The  subscript  i  denotes  an  initial  value.  The  exponent  k  is 
set  by  the  initial  and  final  measured  lengths  and  diameters  of  the 
sample.  A  typical  value  for  k  is  0.85.  The  assumption  of  Eq.  (2) 
is  essentially  arbitrary,  and  other  relations  might  be  suggested, 
e.g., 

L(p,20)-L(p„20)  d{p,20)-d(p„20) 

KP/,20)-Kp„20)  d(p,,20)-d(p„20) 

where  the  subscript  /  denotes  a  final  value.  The  maximum  differ¬ 
ence  in  calculated  sample  lengths  given  by  Eqs.  (2)  and  (3),  for 
any  sample  length,  is  about  0.013  mm.  This  corresponds  to  an 
error  in  density  of  less  than  0.04%.  This  is  considered  negligible 
in  our  experiments. 

The  initial  gap  between  the  sample  and  the  capsule  vanishes 
during  HIP.  Thus  the  sample  diameter  is  the  final  outside 
capsule  diameter  minus  twice  the  final  capsule  wall  thickness,  and 
the  final  sample  length  is  the  final  capsule  length  minus  the  thick¬ 
ness  of  the  end  plugs.  This  assumption  is  supported  by  sample 
examination  after  HIPing.  No  gaps  were  found  after  HIPing.  A  thin 
0.2-micron  reaction  layer  consisting  of  oxides  of  iron,  aluminum, 
chromium,  and  nickel  was  detected  between  the  alumina  sample 
and  the  stainless  steel  capsule.  An  error  in  density  of  less  than 
0.05%  results  from  ignoring  this  layer. 

Putting  all  these  factors  together,  the  program  proceed*  in  the 
following  way.  Pressure,  temperature,  and  shrinkage  are  measured 
for  a  given  time,  and  the  diameter  of  the  capsule  is  calculated  by 
subtracting  the  shrinkage  from  the  original  size.  Then  the  problem 
is  attacked  from  the  other  direction.  A  fractional  density  i*  assumed 
for  the  alumina  sample,  and  the  diameter  and  length  of  the  sample 
at  20°C  are  calculated  using  Eq.  (2)  and  the  mass  of  the  sample. 
The  dimensions  of  the  sample  at  temperature  are  then  calculated 
using  Equation  (1A).  The  volume  of  the  capsule  sidewall  at  20°C 
is  corrected  for  temperature  by  converting  the  linear  expansion  of 
Equation  (IB)  to  a  volume  expansion,  and  the  necessary  wall 
thickness  to  produce  that  volume  is  calculated.  For  the  assumed 
density,  the  outside  diameter  of  the  capsule  will  be  the  sample 
diameter  plus  *wice  the  wall  thickness.  The  program  repeats  this 
process  with  different  values  of  density  until  it  finds  an  outside 
diameter  equal  to  the  actual  diameter  of  the  capsule.  The  output  of 
the  program  is  a  table  consisting  of  five  columns  of  data:  tem¬ 
perature,  pressure,  time,  capsule  shrinkage,  and  fractional  density. 

Ond  difficulty  with  this  approach  is  that  there  are  several  meth¬ 
ods  for  treating  the  sidewall  For  any  of  these  methods,  the  mass 
of  the  sidewall  is  constant  since  it  is  assumed  that  the  sidewall 
does  not  slip  over  the  sample  during  densification,  and  the  mass  is 
given  by 

M  =  itw(D  —  w)Lp  (4) 

where  w  is  the  sidewall  thickness,  D  is  the  outside  diameter  of  the 
capsule,  L  u  the  length  of  sidewall,  and  p  is  the  absolute  density 
of  the  capsule  material.  The  sidewall  mass  can  be  obtained  in  three 
ways:  1)  w  is  assumed  to  be  the  initial  capsule  wall  thickness,  D 
fhe  initial  outside  diameter  of  the  capsule,  and  L  the  initial  capsule 
length  minus  the  initial  length  of  the  end  plugs;  2)  w  is  assumed 
to  be  the  initial  capsule  wall  thickness,  D  the  initial  outside  di¬ 
ameter  of  the  capsule,  and  L  the  initial  sample  length;  and  3)  tv  is 
assumed  to  be  the  final  capsule  wall  thickness,  D  the  final  outside 
diameter  of  the  capsule,  and  L  the  final  samg'e  length.  The  length 
L  used  in  Method  1  is  longer  than  that  used  in  Method  2.  the 
difference  being  the  length  of  the  gaps  between  the  sample  and  the 
end  plugs  (see  Fig.  2).  Method  3  differs  from  the  other  two  since 
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Fig.  3.  Data  for  experiment  No.  2.  (A)  Pressure  (solid  line)  and 
temperature  (dashed  line)  vs  time.  (B)  Relative  density  calculated 
from  temperature  and  capsule  shrinkage  for  alumina  at  1000'C  and 
102  MPa 


Tim*  (min) 

Fig.  4.  Data  for  experiment  No.  4.  (A)  Pressure  (solid  line)  and 
temperature  (dashed  line)  vs  time.  (B)  Relative  density  calculated 
from  temperature  and  capsule  shrinkage.  Density  curve  clearly 
shows  effect  of  changing  pressure  from  34  to  70  MPa  at  330  min. 
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it  uses  final  measurements  instead  of  initial  ones.  Table  II  shows 
the  sidewall  masses  determined  by  these  three  methods.  As  ex¬ 
pected,  the  masses  calculated  by  Method  2  are  consistently  smaller 
than  those  given  by  Method  1,  since  Method  2  assumes  a  slightly 
longer  sidewall.  The  choice  of  sidewall  mass  affects  the  density 
value  calculated  from  the  capsule  shrinkage  data.  Fractional  densi¬ 
ties  can  differ  by  as  much  as  4.6%.  Method  1  can  probably  be 
discounted  since  all  samples,  including  those  that  exhibited  a  small 
density  increase  during  HIPing,  were  bonded  to  the  capsule.  It  is 
likely  that  significant  capsule  deformation  occurred  prior  to  shrink¬ 
age  of  the  alumina  compact  and  removed  the  gap  between  the 
sample  and  the  capsule 

We  have  used  Method  3  to  obtain  the  sidewall  mass  since  it  does 
not  require  any  assumption  regarding  the  nature  of  the  deformation 
of  the  capsule.  In  Table  III,  we  show  the  final  densities  as  calcu¬ 
lated  from  the  probe  data  and  determined  directly  by  the  immersion 
technique.  The  root-mean-squarc  difference  between  the  calcu¬ 
lated  and  measured  densities  is  1.2%.  Among  the  various  possible 
sources  of  error,  the  following  are  probably  the  most  significant. 
First,  there  is  some  uncertainty  m  the  measurements  of  the  initial 
and  final  densities.  These  errors  are  estimated  at  about  0.2%  each 
Second,  there  are  significant  variations  in  the  repotted  values  of  the 
coefficients  of  thermal  expansion,  particularly  for  alumina.  If  the 
coefficient  ot  linear  thermal  expansion  for  alumina  is  in  error  by 
15%,  the  calculated  fractional  density  wili  be  in  error  by  0.3%.  A 


5%  error  in  the  thermal  expansion  of  stainless  steel  would  also 
cause  an  error  of  about  0.3%.  Third,  we  have  assumed  that  the 
sample  remains  cylindrical  during  densification,  but  this  is  not 
quite  correct.  After  HIP,  the  sides  of  the  sample  flare  slightly 
toward  the  ends,  and  the  end  surfaces  become  slightly  domed.  The 
magnitude  of  the  effects  of  this  shape  change  on  tlr  calculated  final 
density  is  unknown.  Finally,  if  there  are  variations  of  density  with 
position  in  the  bisque-fired  sample,  the  dilatometer  may  be  provid¬ 
ing  data  on  a  portion  of  the  sample  which  is  not  typical  of  the  entire 
sample.  If  the  fractional  density  in  the  vicinity  of  the  dilatometer 
probes  differs  from  the  average  for  the  entire  sample  by  0  5%.  the 
calculated  final  density  will  also  be  in  error  by  0.5%.  In  light  of 
these  considerations,  the  calculated  and  measured  densities  appear 
to  be  in  reasonable  agreement. 

Application 

The  dilatometer  and  computer  software  program  enables  volu¬ 
metric  changes  to  be  determined  inside  the  autoclave.  Thus,  we  can 
obtain  quantitative  data  about  tensity  changes  and  phase  changes 
occurring  during  the  HIP  cycle  Reaction  kinetics  data  may  also  be 
obtained  by  measuring  densification  rates  as  a  function  of  density, 
temperature,  and  pressure  and  altering  the  temperature  and' or  pres¬ 
sure  during  the  cycle.  Examples  of  the  use  of  the  dilatometer  are 
given  in  Figs  3  and  4  Figures  3(B)  and  4(B)  show  typical  density  - 
time  plots  for  the  densification  of  alumina  at  1000°C  and 
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Figures  3(A)  and  4(A)  show  the  changes  in  temperature  and  pres¬ 
sure  for  these  two  HIP  cycles  during  pressurization,  heat  up,  and 
HIPing  at  1000°C.  Figure  3(B)  shows  that  significant  densification 
does  not  occur  until  the  temperature  approaches  1000°C,  at  which 
temperature  the  densification  rate  initially  increases  rapidly.  The 
densification  rate  subsequently  falls  with  increasing  fractional  den¬ 
sity  as  expected.  In  experiment  No.  4  (see  Fig.  4)  the  pressure  was 
changed  from  34  MPa  to  70  MPa  during  the  HIP  cycle.  The 
change  in  gradient  of  the  density-time  plot  (Fig.  4(B))  at  a  frac¬ 
tional  density  of  0.72  shows  that  the  HP  densification  of  alu¬ 
mina  is  pressure  dependent.  The  densification  rate  is  proportional 
to  the  square  of  the  pressure.2  The  initial  drop  in  density  shown  in 
Fig.  4(B)  is  not  a  real  effect.  We  believe  it  was  due  to  a  temporary 
movement  of  the  probe  arms  during  initial  gas  pressurization.  It 
was  not  observed  in  any  of  the  other  experiments. 

From  a  practical  standpoint,  considerable  time  and  cost  savings 
are  achieved  using  the  diiatometer.  More  experimental  data  can  be 
obtained  in  one  cycle,  and  the  optimum  conditions  for  HIPing 
powder  compacts  to  full  density  can  be  determined  quickly  by 
altering  the  processing  conditions  during  the  HIP  cycle.  This  ap¬ 
proach  virtually  removes  the  costly  iterative  trial-and-error  method 
commonly  used.  While  the  equipment  has  been  used  mainly  to 
examine  die  densification  of  powder  compacts  it  can  also  be  used 
to  monitor  changes  occurring  in  sintered  ceramics  (e.g.,  densi¬ 
fication,  swelling,  and  phase  changes). 

In  our  work  the  HIP  conditions  have  been  controlled  manually, 
but  extension  to  real-time  process  control  should  be  possible.  A 
primary  obstacle  to  achieving  real-time  control  is  that  a  value  of  k 
is  necessary  in  calculating  density  from  Eq.  (2),  and  k  is  deter¬ 
mined  from  the  final  dimensions  of  the  sample.  Therefore,  it  is 
necessary  to  HIP  at  least  one  sample  to  determine  k  before  real¬ 
time  control  is  possible.  This  appears  to  be  a  reasonable  strategy 
since  our  work  indicates  that  k  does  not  vary  strongly  from  sample 
to  sample. 

To  date,  the  equipment  has  been  used  only  at  relatively  low 
temperatures  (S1200°C)  to  study  the  densification  of  alumina.  The 
stainless  steel  capsules  can  be  used  up  to  1400°C,  but  if  molyb¬ 
denum  or  tantalum  were  used  as  the  capsule  material,  the  equip¬ 
ment  couid  be  used  up  to  1750°C.  The  furnace  would  have  to  be 
modified  in  order  to  operate  above  1800°C. 
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Interface-Reaction-Controlled  Kinetics  in  the 
Hot  Isostatic  Pressing  of  Submicrometer 
Alumina  Powder 

Roger  R.  Wills*  and  John  Kevin  McCoy 

Battelle's  Columbus  Laboratories,  505  King  Avenue,  Columbus,  Ohio  43201-2693 

The  kinetics  of  densification  of  alumina  powder  compacts  at  1000 0  to  1200°C  were 
studied  by  using  a  dilatometer  to  monitor  volumetric  changes  during  hot  isostatic 
pressing.  The  densification  rate  is  proportional  to  the  square  of  the  applied 
pressure,  and  densification  does  not  occur  below  a  threshold  pressure.  Results  are 
interpreted  in  terms  of  an  interface -controlled  reaction  mechanism  involving  the 
movement  of  grain-boundary  dislocations. 
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A  n  empirical  approach  is  generally  used 
**to  determine  the  optimum  conditions 
for  hot  isostatic  pressing  of  ceramic  powder 
compacts  to  full  density.  We  have 
attempted  to  remove  this  empiricism  by 
examining  mechanisms  responsible  for 
densification  during  hot  isostatic  pressing. 
Furthermore,  by  comparing  the  predictions 
of  various  constitutive  equations  with  ex¬ 
perimental  data,  we  have  attempted  to 
show  that  optimum  hot  isostatic  pressing 
conditions  can  be  forecast  if  certain  infor¬ 
mation  about  the  powder  is  known.  Alumi¬ 
num  oxide  was  selected  for  study  since 
there  is  a  wealth  of  data  cm  this  ceramic. 
The  present  communication  presents  evi- 
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dence  to  indicate  that,  under  hot  isostatic 
pressing  conditions,  the  densification  of 
undoped  alumina  is  controlled  by  an  inter¬ 
face  reaction  involving  the  movement  of 
grain-boundary  dislocations.  We  believe 
that  this  is  the  first  time  that  this  mechanism 
has  been  observed  to  be  active  in  the  densi¬ 
fication  of  a  ceramic  powder  compact. 

Previous  research1*2  suggested  that 
the  mechanism  of  densification  of  sub- 
micrometer  alumina  powder  was  probably 
grain-boundary  diffusion.  After  the  initial 
stage  of  particle  rearrangement,  densi¬ 
fication  proceeds  by  grain-boundary  dif¬ 
fusion,  causing  the  average  number  of 
contacts  per  particle  to  increase.  Arzt 
et  al. 3  expressed  this  stage  of  densification 
mathematically  as’ 


dD/dt= 


l28DbaZP* 

rMD)kT 


(1) 


where  S  is  the  effective  grain-boundary 
thickness,  Db  the  grain-boundary  diffusion 
coefficient,  ft  the  atomic  volume,  P*  the 
driving  pressure  for  densification,  k  is 
Boltzmann’s  constant,  T  temperature,  Z 


the  particle  coordination  number,  h(D)  a 
geometric  factor  depending  on  density,  and 
r  the  particle  radius.  The  equation  predicts 
that  the  rate  of  densification  is  proportional 
to  the  driving  pressure  P  * .  This  is  related  to 
the  externally  applied  pressure  as  follows: 

P*=^~-+P-P.  (2) 

D 

where  P  is  the  externally  applied  pressure, 
F,  the  effective  contact  pressure  resulting 
from  surface  tension,  P,  the  internal  pore 
pressure,  D  the  fractional  density  of  the 
compact,  and  a  the  average  neck  area. 
Since  powder  compacts  are  omgassed 
before  hot  isostatic  pressing,  P,  is  =0. 
Furthermore,  P,  can  be  considered  to  be 
negligible  in  comparison  with  the  first  term 
in  Eq.  (2),  particularly  in  view  of  the  rela¬ 
tively  high  pressures  used  in  pressing.  At 
constant  density,  die  net  driving  pressure  is 
consequently  proportional  to  the  applied 
pressure.  Thus,  by  varying  the  pressure 
during  the  hot  isostatic  pressing  cycle,  the 
dependency  of  the  rate  of  densification  on 
the  external  pressure  can  be  determined. 
Figure  1  shows  the  effect  of  changing  the 
pressure  from  34  to  70  MPa  on  the  densi¬ 
fication  of  a  powder  compact  encapsulated 
in  a  stainless-steel  capsule.  Assuming  that 
the  rate  of  densification  is  proportional  to 
Pn,  m  is  =2.1.  However,  if  grain-bound¬ 
ary  diffusion  is  controlling  the  densification 
rate,  the  exponent  should  be  equal  to  unity 
(see  Eq.  (1)).  One  possible  explanation  for 
the  value  of  the  pressure  exponent  is  that  a 
combination  of  mechanisms  is  operative 
during  hot  isostatic  pressing.  For  example, 
a  dislocation-glide  mechanism  gives  a  pres¬ 
sure  exponent  of  3/  and  consequently,  if  it 
operated  in  conjunction  with  grain¬ 
boundary  diffusion,  a  pressure  exponent  of 
2.1  might  result.  This  explanation  can  be 
discounted  since  the  climb  or  glide  of  lat¬ 
tice  dislocations  is  unlikely  to  occur  in  the 
temperature  range  1000°  to  1200°C.  Anoth- 
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Fig.  1.  Effect  of  changing  applied  pressure  from  34  to  70  MPa  on 
denstfication  of  alumina  at  fractional  density  of  0.725  (1000°C). 


Fig.  2  Densification  rate  vs  applied  pressure  for  fractional  densities 
indicated. 


Table  1.  Threshold  Stress  Values  at  Different  Temperatures  and  Fractional  Densities 


D 

Temperature  (*C) 

n 

P.  (MPa) 

P»  (MPa) 

0.71 

1050 

6.2478 

24.435 

80.66 

0.78 

1050 

2.9318 

8.564 

61.33 

0.71 

1000 

6.2478 

38.645 

132.35 

0.72 

1000 

5.4152 

36.403 

101.38 

er  possibility  is  that  basal-plane  slip  is  re¬ 
sponsible  for  densification.  This  mech¬ 
anism  may  contribute,  but  it  could  not  be 
the  only  factor  since  it  does  not  enable  mass 
transport  to  occur  in  all  directions. 

One  mechanism  that  would  explain 
die  results  is  interface-reaction-controlled 
grain-boundary  diffusion.  It  is  considered 
to  be  operative  in  many  creep  experiments 
at  low  temperatures  and  small  grain  sizes 
where  it  has  been  suggested  that  the  grain 
boundaries  do  not  act  as  perfect  sources  or 
sinks  of  point  defects/  and  the  rate  of  the 
interface  reaction  that  creates  and  anni¬ 
hilates  vacancies  controls  creep.  Ashby’ 
suggested  that  vacancy  creation  and 
annihilation  occur  by  the  climb  of  grain¬ 
boundary  dislocations.  Burton*  assumed 
that  spiral  dislocation  climb  sources  main¬ 
tained  the  dislocation  density  as  creep  pro¬ 
ceeded  and  dislocations  migrated  along  the 
grain  boundary. 

Adapting  the  treatment  of  Arzt  et  al. 7 
to  densification,  one  finds  that  the  densi¬ 
fication  rate  is  proportional  to  the  square  of 
the  driving  pressure  and  that  the  expression 
for  the  driving  pressure  is  modified  by  the 
addition  of  a  threshold  term,  P Below 
this  threshold  pressure,  there  is  insufficient 
sues;  to  cause  dislocation  motion.  Modi¬ 
fying  Eq.  (1)  for  the  hot  isostatic  pressing 
of  powder  compacts  gives: 


P*  = 


Airr3P 

aZD 


+P-P* 


(3) 


The  existence  of  a  threshold  stress  can  be 
shown  from  the  experimental  data.  At  con¬ 
stant  density,  temperature,  and  grain  size 


dD/dt=QP*} 


(4) 


where  Q  is  a  constant.  Equation  (3)  indi¬ 
cates  that  densification  can  only  occur  by 
this  interface  mechanism  if  P*  is  >0  or 
( Airr*P/aZD  +/’,)>/’, *.  When  no  external 
pressure  is  applied,  densification  is  not  ob¬ 
served.  Thus,  P,h  must  be  greater  than  or 
equal  to  P„  otherwise  some  densification 
would  be  observed  in  this  sintering  regime. 
Now 

dDldt=QP*3=^y~£-+P-P^  (5) 
Patting  n  =4nr3/oZD ,  and  P,-Pu,-P„ 

(< £?/dtya=nQ'a(p (6) 

If  the  square  root  of  the  densification  rate 
is  plotted  against  applied  pressure,  the 
graph  should  intercept  the  positive  x-axis. 
Figure  2  shows  the  graphs  for  four 
conditions:  (1)  D=0.71,  T=1050°C, 
(2)  D=0.78,  r=1050°C,  (3)  D=0.71, 
r=1000°C,  and  D=0.72.  r=1000°C.  At 
zero  densification  rate,  the  plots  intercept 
the  pressure  axis  as  predicted.  Since  the 
intercept  value  is  equal  to  (Prt-F,)/n ,  ex¬ 
act  values  for  P*  cannot  be  obtained  from 
these  data  without  knowledge  of  P,  for 
these  four  conditions.  These  were  calcu¬ 
lated  using  the  equations  derived  by  Aizt 
and  co-workers.3'*  Table  I  shows  the  values 
for  the  threshold  stress.  As  expected,  these 
values  decrease  with  both  increasing  frac¬ 
tional  density  and  increasing  temperature. 

Interface-controlled  diffusional  creep 
has  been  found  in  AI2Oi/  U02,'°  Fe203, 
and  MgO. 12  For  example.  Cannon  et  ol.'> 
reported  a  threshold  stresc  of  13.8  MPa  for 


1.2-/im  grain  size  A1203  doped  with 
0.25%  MgO.  The  strain  rate  also  varied 
with  the  square  of  the  stress  for  both  pure 
and  MgO-doped  AI203  at  strain  rates  below 
10"’  s'1.  These  data  suggest  that  the 
mechanism  of  densification  may  also  be 
interface-controlled  in  the  latter  stage  of  the 
densification  (D>0.90)  of  alumina  and 
therefore  could  be  the  rate-controlling 
mechanism  over  most  of  the  densification 
range  (D=0.65  to  1.0). 
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Fig  9— Metallography  section  of  impacted  semi-infinite  plate  of  hard¬ 
ness  150  HV.  Impaci  velocity  325  ms'1  Magnification  6  8  times. 


formation  of  a  dead  zone  below  the  compressed  region.  In 
comparing  Figures  1  and  9,  it  is  apparent  that  a  prerequisite 
for  shear  band  initiation  was  some  degree  of  acceleration  of 
the  section  of  plate  between  the  advancing  projectile  and  the 
distal  boundary  of  the  target,  consistent  with  the  slip  line 
field  model  proposed  by  Packman  and  Finnegan2  which 
employed  a  minimum  plastic  strain  rate  criterion  for  shear 
band  formation.  The  shear  stress  concentrating  effect  of  the 
projectile  geometry  was  therefore  not  in  itself  sufficient  for 
shear  band  initiation. 

This  work  was  carried  out  at  the  Materials  Research 
Laboratories,  Maribymong,  Victoria,  Australia. 
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Discussion  of  “Practical  Applications 
of  Hot-Isostatic  Pressing  Diagrams: 
Four  Case  Studies”* 

J.  KEVIN  McCOY 

A  recent  paper  by  Arzt,  Ashby,  and  Easterling1  describes 
a  theory  of  densification  by  hot-isostatic  pressing.  The  pur¬ 
pose  of  this  note  is  to  point  out  an  error  that  was  made  early 
in  the  mathematical  development  which  affected  some  of  the 
subsequently  obtained  results. 

Arzt,  Ashby,  and  Easterling  use  two  dis'met  geometrical 
descriptions  of  densification.  Since  these  two  descriptions 
are  mathematically  equivalent,  Arzt  et  al.  switch  from  one 
to  the  other  as  best  suits  the  aims  of  clarity  and  mathematical 
convenience.  However,  it  is  necessary  to  keep  firmly  in 
mind  which  description  is  currently  in  use. 

The  first  description  I  will  call  the  “growing-sphere 
model".  This  is  a  relatively  new  approach,  having  been 
pioneered  by  Arzt.2  The  growing-sphere  model  is  illustrated 
in  Figure  1:  a  group  of  spherical  particles  is  assumed  to 
grow  around  fixed  centers.  As  the  particles  grow,  they  will 
impinge  on  one  another.  Arzt2  has  treated  this  effect  in  two 
ways;  following  Arzt,  Ashby,  and  Easterling,1 1  will  discuss 
only  one  of  these.  In  this  model,  it  is  assumed  that  matter  is 
forced  out  of  the  contact  region  to  form  a  localized  neck, 
while  away  from  necks,  the  particle  radius  and  spherical 
geometry  are  unaltered.  Therefore,  if  the  particles  grow 
from  a  radius  R  to  a  radius  /?’,  the  volume  of  each  particle 
increases  from  4ttR 3/3  to  4-rrR  ’3/3.  Since  the  particles  grow 
around  fixed  centers,  the  relative  density  increases  from  its 
initial  value  D0  to  D,  where 


This  is  equivalent  to  Eq.  [  1  ]  of  Arzt,  Ashby,  and  Easterling. 

The  second  description  I  will  call  the  “colliding-sphere 
model”.  This  is  the  standard  description  of  densification. 
used,  for  example,  by  Wilkinson  and  Ashby.3  The  colliding- 
sphere  model  is  illustrated  in  Figure  2.  In  this  treatment,  the 
spherical  panicles  are  assumed  to  have  fixed  volume,  but 
the  particle  centers  are  assumed  to  move  so  that  each  panicle 
approaches  every’  other  particle.  One  might  say  that  the 
particle  volumes  remain  constant  while  space  shrinl  s.  This 
wording  points  out  the  mathematical  equivalence  of  the 
growing-sphere  and  colliding-sphere  models,  they  differ 
only  by  a  dilatation  of  space. 

Using  the  colliding-sphere  model,  suppose  that  densi- 
fication  reduces  the  c»  nter-to-center  distance  of  two  initially 
tangent  particles  from  2 R  to  2y  while  increasing  the  relative 
density  from  D0  to  D.  Since  the  colliding-sphere  model 
assumes  a  uniform  contraction  of  the  group  of  particles. 


*E  ARZT.  M  F  ASHBY,  and  K  E  EASTERLING  Metall  Trans  A. 
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DENSIFICATION  BY  INTERFACE-REACTION  CONTROLLED  GRAIN-BOUNDARY  DIFFUSION 


J.  KEVIN  McCOY  and  ROGER  R.  WILLS 
Battelle,  505  King  Avenue,  Columbus,  Ohio  43201-2693,  USA 


Abstract — The  denslflcatlon  of  a  fine-grained,  high-purity  aluminum 
oxide  powder  under  hot  Isostatic  pressing  (HIP)  has  been  found  to  occur 
by  Interface-reaction  controlled  grain-boundary  diffusion.  We  discuss 
geometries  and  dislocation  mechanics  for  this  process  for  both  the  ini¬ 
tial  and  final  stages  of  denslflcatlon  and  develop  constitutive  equa¬ 
tions  for  denslflcatlon  rate  as  a  function  of  density,  materials  con¬ 
stants,  and  experimental  parameters.  The  model  Is  used  to  explain  the 
results  of  several  HIP  experiments  at  pressures  of  34  to  1C2  MPa  and 
temperatures  of  1273  to  1423  K.  Sources  of  variation  from  sample  to 
sample  are  discussed.  An  analysis  Is  made  of  the  sensitivity  of  the 
model  to  Its  adjustable  parameters.  Alternative  explanations  for  the 
experimental  data  are  discussed  and  found  to  be  Inadequate. 


1.  INTRODUCTION 

A  number  of  mechanisms  have  been  proposed  to  describe  denslfl¬ 
catlon  during  sintering  and  hot  Isostatic  pressing  (HIP).  Most  of  the 
theoretical  treatments  Invoke  either  a  simple  diffusion  or  dislocation 
mechanism  for  material  transport.  When  we  analyzed  densiflcation  data 
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for  aluminum  oxide  from  a  series  of  HIP  experiments  In  which  both  tem¬ 
perature  and  pressure  were  monitored*  we  found  that  the  data  could  not 
be  adequately  explained  In  terms  of  either  theory*  but  that  they  could 
be  explained  In  terms  of  Interface- react Ion  controlled  grain-boundary 
diffusion  Cl],  Both  dislocation  and  diffusion  elements  are  contained  In 
this  mechanism  in  that  mass  transport  occurs  by  diffusion  but  Is  limited 
by  the  mobility  of  grain-boundary  dislocations.  The  concept  of 
Interface-reaction  controlled  grain-boundary  diffusion  Is  not  particu¬ 
larly  new,  having  been  proposed  by  Burton  [2]  In  1972  and  giver,  a 
thorough  development  In  a  recent  review  by  Arzt*  Ashby,  and  Verrall  [3], 
but  It  has  apparently  been  applied  only  to  creep  or  deformation. 

In  this  paper,  we  apply  the  theory  of  Interf ace- react Ion  con¬ 
trolled  grain-boundary  diffusion  to  describe  denslflcatlon  of  a  porous 
solid.  Following  standard  practice  [4],  we  divide  denslflcatlon  Into 
two  stages.  The  initial  stage  Is  described  as  a  random  packing  of 
spherical  particles;  this  model  Is  applied  ax  fractional  densities  of 
0.64  to  0.9.  The  final  stage  Is  described  by  an  Isolated  spherical  pore 
centered  In  a  spherical  shell;  this  model  Is  used  for  fractional  densi¬ 
ties  greater  than  0.9.  We  develop  appropriate  geometries  for  both 
stages  of  denslflcatlon  and  then  briefly  review  the  islocation  kinetics 
associated  with  the  mechanism.  The  model  is  expressed  In  a  form  suit¬ 
able  for  fitting  to  experimental  data  and  applied  to  explain  the  denslf- 
Icatlon  kinetics  of  a  high-purity  aluminum  oxide  for  a  range  of  tempera¬ 
tures,  pressures,  and  specimen  porosities. 
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2.  INITIAL  STAGE 


In  this  section,  we  develop  a  model  for  the  Initial  stage  of 
denslflcatlon.  We  begin  by  describing  the  geometry  of  a  neck  between 
two  particles  and  use  the  results  to  calculate  the  driving  force  for 
denslflcatlon.  Finally,  the  drlvlnq  force  and  the  mechanism  of 
Interface-reaction  controlled  grain-boundary  diffusion  are  used  to  cal¬ 
culate  the  denslflcatlon  rate. 

The  most  Important  geometrical  quantities  for  the  description 
of  an  Interparticle  neck  are  Illustrated  In  Fig.  1:  R  Is  the  particle 
radius,  2y  Is  the  particle  center-to-center  distance,  x  Is  the  neck 
radius,  and  w  Is  the  axial  radius  of  curvature  of  the  neck.  These  quan¬ 
tities  reflect  a  standard  model  of  neck  geometry:  twc  spherical  parti¬ 
cles  impinge  upon  each  other,  and  matter  Is  transferred  from  the  contact 
area  to  a  localized  neck  while  the  remainder  of  the  surface  of  the 
spheres  remains  unaltered.  The  surface  of  the  neck  Is  assumed  to  be  a 
portion  of  a  circular  torus  that  Is  tangent  to  both  spheres.  The  minor 
radius  of  the  torus  Is  w  and  the  major  radius  is  x+w.  For  given  values 
of  R  and  y,  x  and  w  are  specified  by  the  solution  to  two  equations.  The 
first  equation  follows  from  the  Pythagorean  theorem: 

y2+(x+w)2  =  (R+w)2  (1) 

or 


. . 

w  “  2(R-x)  • 

The  second  equation  for  x  and  w  specifies  that  the  volume  of  matter 
removed  from  the  spheres  must  equal  the  volume  of  matter  deposited  in 


(2) 


4 


the  neck.  It  Is  a  straightforward  exercise  In  calculus  to  show  that 
;h1s  equation  Is 

|(R-y)2(2R+y)  =  (3) 

R2wcos20s1n0  -  ws1n0(R2-y2) 

+  w^sln^©  -  2(x+w)w2(0-cos6s1n0) 

where 


0  =  arctan[y/(x+w)3.  (4) 

Due  to  the  complexity  of  equation  (3),  Arzt  [53  has  proposed  the  approx¬ 
imation 

x2  *  n(R-y)R  (5) 

where  n  Is  a  numerical  constant.  Arzt  suggested  the  value  n  =  II,  but 
we  have  chosen  to  use  n  =  11.4  since  this  value  gives  a  slightly  better 
fit.  We  have  used  equations  (2)  and  (5)  to  provide  a  description  of 
torus  geometry.  As  evidence  of  the  accuracy  of  equation  (5),  we  give  In 
Fig.  2  a  comparison  of  the  values  of  x  and  w  obtained  by  using  the 
rigorous  treatment  of  equation  (3)  and  the  approximation  of  equation 
(5).  The  agreement  between  the  two  treatments  Is  probably  better  than 
that  between  either  treatment  and  physical  reality. 

With  this  description  of  neck  geometry,  we  can  find  p*,  the 
effective  contact  pressure  between  the  two  particles  for  the  Initial 
stage,  which  produces  a  climb  force  on  grain-boundary  dislocations  and 
thus  acts  as  the  driving  force  for  denslflcatlon.  Following  Molerus  [6] 
and  Arzt,  Ashby,  and  Easterling  [73,  we  write  the  effective  contact 


pressure  as 
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5*  =  Asfiji 


rfDVth 


(6) 


where  p  Is  the  externally  applied  pressure#  a  Is  the  average  neck  area# 

Z  Is  the  coordination  number#  D  Is  the  fractional  density#  p  Is  the 
contact  pressure  due  to  surface  tension#  and  p.^  Is  a  threshold  pressure 
for  the  motion  of  grain-boundary  dislocations.  In  our  experiments#  the 
sample  was  thoroughly  outgassed  before  encapsulation#  and  the  capsules 
were  welded  in  vacuum#  so  the  effect  of  gas  trapped  In  the  pores  Is  not 
Included  In  equation  (6). 

To  evaluate  the  effects  of  coordination  number  and  contact 
area  In  equation  (6)»  It  Is  mathematically  convenient  to  use  the 
approach  of  Arzt  [5].  Arzt  describes  denslflcatlon  not  as  a  packing  o4- 
spheres  that  approach  each  other  but  as  a  packing  of  stationary  spheres 
that  grow  fictitiously.  If  the  spheres  grow  from  a  radius  R  to  a  radius 
R'»  the  Increase  In  the  volume  of  the  spheres  will  cause  the  fractional 
density  to  Increase  from  Its  initial  value  to  a  new  value  D,  where 


D/DQ  =  (R»/R)3  (7) 

As  the  spheres  grow#  two  types  of  geometrical  changes  occur. 
First#  the  Initial  contact  points  grow  Into  necks  as  the  growing  spheres 
Impinge  on  each  other.  We  have  assumed  that  all  displaced  material 
remains  In  the  necks#  while#  away  from  the  necks#  the  particles  remain 
spherical  with  radius  R’.  Second#  as  sphere  growth  continues#  the 
spheres  also  begin  to  form  contacts  with  additional  spheres  that  they 
did  not  touch  in  the  original  packing.  Using  radial  distribution  func¬ 
tions  for  a  random  dense  packing  developed  by  Scott  C8]  and  Mason  [9]# 
Arzt  developed  the  following  approximate  equation  for  the  coordination 


number 


Z  =  Z0+c(RVR-l)  (8) 

where  ZQ  =  7.3  and  c  =  15.5.  Using  this  radial  distribution  function 
and  equation  (5)  for  the  area  of  a  neck,  the  average  neck  area  Is  [53 

a  =  n(Zo(R»-R)R+c(R’-R)2/2)R/ZR*.  (9) 

Equation  (9)  Includes  a  correction  so  that  It  Is  applicable  to  spheres 
of  radius  R  rather  than  R’.  Using  equations  (8)  and  (9),  the  effect  of 
external  pressure  In  equation  (6)  may  be  evaluated. 

Following  Arzt,  Ashby,  and  Easterling  [73,  the  contact  pres¬ 
sure  due  to  surface  tension  Is 

(10) 

where  y  Is  the  surface  tension.  The  average  neck  area  Is  known  from 
equation  (9),  so  the  average  neck  radius  Is 

x  =  Va/n  (11) 

while  the  axial  radius  of  curvature  of  the  neck  Is  given  by  equation 

(2). 

The  threshold  pressure  for  motion  of  a  grain-boundary  disloca¬ 
tion  Is  uncertain,  since  It  depends  on  the  Burgers  vector  of  the  dislo¬ 
cations  and  the  type  of  dislocation  source.  We  have  chosen  to  use  the 
approach  of  Arzt,  Ashby  and  Verrall  [33.  Following  their  treatment,  the 
threshold  stress  for  grain-boundary  dislocation  motion  Is 


(12) 


where  a  =  0.2  Is  a  constant  that  reflects  the  relative  amount  by  which 
o 
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the  length  of  a  dislocation  fluctuates  as  It  moves. 

The  presence  of  an  effective  contact  pressure  will  cause  den- 
si  flcatlon  and  produce  changes  In  the  neck  geometry.  As  densl flcatlon 
proceeds#  the  center-to-center  distance  of  two  Initially  tangent  parti¬ 
cles  will  change  with  rate  d2y/dt.  As  discussed  below  In  the  section  on 
dislocation  mechanics#  a  change  In  y  results  from  the  motion  of  grain¬ 
boundary  dislocations;  the  passage  of  one  dislocation  across  the  entire 
neck  reduces  the  center-to-center  distance  2y  by  b^  where  b  Is  the  com- 
ponent  of  the  Burgers  vector  normal  to  the  grain  boundary.  For  the  case 
In  which  several  dislocations  are  moving  together#  the  number  cf  dislo¬ 
cations  that  will  pass  a  given  point  on  the  grain  boundary  In  unit  time 
Is  vp.  Thus  the  rate  of  change  of  the  center-to-center  distance  of  two 
particles  Is 


*  -'"V  <»> 

The  final  step  of  the  model  Is  to  obtain  the  densl flcatlon 
rate.  If  It  Is  assumed  that  the  sample  starts  as  a  packing  of  tangent 
spherical  particles  and  that  the  particles  approach  each  other  as  den- 
slflcatlon  proceeds#  It  Is  clear  from  geometry  that  [103 


D  =  DQ(R/y)3  (14) 

where  Dq  is  the  Initial  fractional  density  and  D  Is  the  current  frac¬ 
tional  density.  By  differentiating  equation  (14)#  one  obtains 


1/3 


An  expression  for  dD/dt  In  termt  of  experimental  variables  Is  given  In 
equation  (27). 
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3.  FINAL  STAGE 


As  a  material  approaches  full  density,  It  Is  no  longer 
appropriate  to  describe  Interparticle  contacts  In  terms  of  distinct 
necks,  and  It  becomes  appropriate  to  describe  the  porosity  In  terms  of 
Isolated  pores.  In  this  section,  we  develop  a  model  for  this  stage  of 
denslflcatlon.  We  first  describe  the  geometry  of  the  grains  and  pores, 
then  determine  the  driving  force  for  denslflcatlon,  and  finally  obtain 
the  denslflcatlon  rate. 

The  geometry  of  the  pores  and  grains  Is  similar  to  that  used 
by  Arzt,  Ashby,  and  Easterling.  It  Is  assumed  that  the  grains  take  the 
shape  of  truncated  octahedrons  with  spherical  pores  of  radius  r  located 
at  each  vertex.  Each  truncated  octahedron  has  24  vertices,  and  four 
grains  meet  at  each  vertex  so  there  are,  on  the  average,  six  pores  per 
grain.  The  volume  of  each  grain  Is  assumed  to  be  4nR^/3,  as  In  the  Ini¬ 
tial  stage.  Since  there  are  six  pores  per  grain,  the  density,  pore 
radius  and  equivalent  spherical  grain  radius  are  related  by 


1-D  = 


(16) 


rW 

Following  Coble  [11],  we  may  simplify  the  geometry  for  the  purposes  of 
calculating  denslflcatlon  rate  by  taking  a  spherical  shell  of  material 
around  each  pore.  We  assign  a  volume  of  matter  2nR^/9  (1/6  of  a  grain) 


to  each  pore,  so  the  shell  radius  Rg  Is  given  by 


3  ,  B^+r3 

s  6  +r 


(17) 
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The  effective  pressure  or  driving  force  for  denslflcatlon  may 
be  obtained  by  balancing  the  forces  Imposed  on  the  shell  against  the 
stresses  In  the  shell.  Suppose  that  we  make  an  Imaginary  cut  through 
the  center  of  the  shell  on  the  xy  plane.  The  total  force  on  one  hemi¬ 
sphere  In  the  z  direction  Is 

jiR2P  +  2nry  -  Jt(R2-r2)c00  =  0.  (18) 

5  S  wo 

The  first  term  Is  due  to  the  externally  applied  pressure*  the  second  Is 
from  the  surface  tension  of  the  pore*  and  the  third  results  from  the 
tangential  stress  within  the  sphere  which  balances  the  forces  applied  to 
the  surface  of  the  sphere.  Note  that  Is  treated  as  a  constant  in 
equation  (18).  This  Is  justified  If  we  assume  that  the  motion  of  dislo¬ 
cations  Is  at  steady  state.  As  discussed  below*  the  density  and  speed 
of  the  dislocations  are  both  proportional  to  the  driving  force.  But  at 
steady  state,  vp  Is  Independent  of  position,  so  cr0e  Is  also  Independent 
of  position. 

In  addition  to  the  tangential  stress,  the  driving  force  for 
denslflcatlon  also  Involves  a  threshold  term*  just  as  for  the  Initial 
stage.  For  hot  Isostatic  pressing*  this  Is  [3] 


\/2a  GbK 
_  _  _ Q-l.fr 

pth  ”  2R 

and  the  net  driving  force  Is 


(19) 


P*  ’  "ee-Pth-  <2° 

Differentiating  equation  (16)  with  respect  to  time,  we  obtain 

dfi  _  =2P_0=BI  dn 

dt  "  r  dt  * 


(21) 


Furthermore* 


dn _ L_  di 


where  V  is  the  volume  of  a  pore  and  -dV/dt  is  the  rate  at  which  matter 


Is  removed  from  the  grain  boundaries  and  deposited  in  the  porej 


dV/dt  =  -vpbnn(R^-r2)nb,  (23) 

that  Is*  the  product,  of  the  dislocation  flux*  the  component  of  the 
Burgers  vector  normal  to  the  grain  boundary*  the  area  of  one  planar 
grain  boundary  Intersecting  the  pore,  and  the  effective  number  of  boun¬ 
daries.  An  expression  for  dD/dt  In  terms  of  experimental  variables  Is 
given  In  equation  (30). 


4.  DISLOCATION  MECHANICS 

In  both  the  Initial  and  final  stages,  the  applied  pressure 
results  In  a  climb  force  on  the  grain-boundary  dislocations.  In  this 
section,  we  discuss  dislocation  effects  that  are  common  to  both  stages. 

Following  Arzt,  Ashby*  and  Verrall  [33*  the  rate  at  which  an 
Individual  grain  boundary  dislocation  climbs  Is 

v  =  MF  =  Mp*bn.  (24) 

Here  F  Is  the  force  per  unit  length  on  the  dislocation  and  M  Is  the 
dislocation  mobility.  From  elementary  dislocation  theory,  it  follows 
that  the  force  on  the  dislocation  Is  the  product  of  the  driving  force  p* 
anti  the  component  of  the  Burgers  vector  normal  to  the  grain  boundary  bn. 
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Applying  the  results  of  Cottrell  and  Jaswon  [12]  and  Cottrell 
[13] ,  Arzt,  Ashby,  and  Verrall  obtained  the  following  expression  for 
grain-boundary  dislocation  mobility  limited  by  solute  drag: 


PkTb^o 


In  this  equation,  Cp  Is  the  concentration  of  solute  In  the  lattice,  and 
the  concentration  In  the  solute  atmosphere  of  the  dislocation  Is 
Increased  by  a  factor  of  p  to  pCQ.  Dg  Is  the  dlffuslvlty  of  the  solute 
In  the  lattice,  Q  Is  the  effective  atomic  volume  of  the  solute,  kT  has 
Its  usual  meaning,  and  b^  Is  the  Burgers  vector  of  the  grain-boundary 
dislocation. 

To  complete  the  description  of  a  grain-boundary  with  disloca¬ 
tions  we  need  only  an  expression  for  the  density  of  dislocations.  This 
problem  has  been  studied  by  Burton  [2]  and  by  Arzt,  Ashby,  and  Verrall 
[3].  Although  different  approaches  were  used,  their  results  are  In  sub¬ 
stantial  agreement.  The  result  of  Arzt  et  al .  Is 


C,p* 


where  C^  Is  a  constant  (about  0.5),  G  Is  the  shear  modulus,  and  bb  Is 
the  length  of  the  Burgers  vector. 


5.  EXPERIMENTAL  RESULTS 


The  development  of  the  theory  of  denslflcatlon  presented  above 
was  prompted  by  failure  of  standard  theories  to  explain  the  data 
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reported  below.  The  experimental  technique  used  In  obtaining  the  data 
has  been  described  In  detail  elsewhere  [14],  so  only  a  brief  summary 
will  be  given  here.  Cylindrical  samples  of  Reynolds  RC-HP-DBM  alumina 
were  encased  In  close-fitting  AISI  316  stainless  steel  capsules.  The 
samples  were  subjected  to  HIP  at  temperatures  ranging  from  1273  K  to 
1423  K  and  at  pressures  ranging  from  34  MPa  to  102  MPa.  Throughout  den- 
slflcatlon,  continuous  records  of  temperature  and  pressure  were  kept. 

In  addition,  a  two-probe  dllatometer  was  used  to  produce  a  continuous 
record  of  the  diameter  of  the  sample.  The  records  of  temperature  and 
diameter,  along  with  a  description  of  the  changes  In  geometry  of  the 
sample,  were  used  to  calculate  a  record  of  density  as  a  function  of  time 
during  HIP.  This  procedure  produces  far  more  Information  from  each 
experiment  than  does  conventional  HIP  practice,  which  provides  only  a 
final  density.  In  addition  to  a  providing  a  much  larger  body  of  data, 
the  dllatometer  technique  provides  the  additional  advantage  that  all  the 
data  obtained  from  a  given  experiment  can  be  compared  to  theory  without 
consideration  of  sample-to-sample  or  run-to-run  ''arlatlons.  To  further 
Increase  the  amount  of  Information  produced  during  each  experiment,  the 
temperature  and/or  pressure  were  often  changed  during  an  experiment. 

This  allowed  us  to  determine  activation  energies  and  pressure  exponents 
from  a  single  experiment. 

For  the  purposes  of  fitting  the  theory  to  the  experimental 
data.  It  is  convenient  to  collect  the  variables  that  are  essentially 
constant  throughout  the  experiment.  If  this  Is  done,  we  may  combine 
equations  14,  16,  and  24  to  26  and  write  for  the  Initial  stage 


t  — j\m  VIIAJI.OA  s^ir  dt'  n*  re^J'tA  ’VR7V*1  -TiJK.  TV*  Tin  /  Ul.  /tXJtX  T^M  TJLntUS 
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1/3  p*2e>:p<-E5/12r) 

GT 


where 


(27) 


f 

2R  [  bb  J  bbPkC0 

and  the  solute  dlffuslvlty  Is  assumed  tc  be  given  by 


(28) 


Ds  =  D°exp(-Eg/RT)  (29) 

Similarly#  we  may  combine  equations  21  to  26  and  write  for  the  final 
stage 

£.  tf]1/V/3(:-(l-D)2/3)p.2^-!-^nbf  (30) 

using  the  same  definition  for  f.  The  following  data  were  used  with 

2 

equations  (27)  and  (30).  The  surface  energy  was  taken  to  be  r  =  1  J/nr 

[73.  The  activation  energy  for  solute  diffusion  Eg  was  determined 

experimentally  by  changing  the  temperature  during  the  course  of  the  run. 

A  value  of  E  =  290000  J/mol  was  measured.  Values  of  the  shear  modulus 
s 

G  were  obtained  from  the  equation 

G  =  12.382X109(1838.5-T)1/3  Pa,  (31) 

where  T  is  the  temperature  In  kelvlns,  which  provides  a  good  fit  to  the 
data  of  Ryshkewltch  [153.  The  value  of  n^,  the  effective  number  of 
grain  boundaries  meeting  at  a  pore,  was  derived  from  the  assumption  that 
the  grain  shape  Is  a  truncated  octahedron,  Jn  this  geometry,  four  regu¬ 
lar  hexagonal  faces  and  two  square  faces  meet  at  each  vertex,  so  th*» 
total  of  the  angles  of  all  the  faces  meeting  at  the  vertex  is  660*. 

Since  a  single  planar  grain  boundary  corresponds  to  an  angle  of  360*,  n^ 
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=  11/6. 

With  these  data#  It  Is  possible  to  use  the  temperature  and 
pressure  records  to  Integrate  the  kinetic  equations  for  densl flcatlon  to 
obtain  a  theoretical  record  of  density  as  a  function  of  time.  Results 
of  such  calculations  are  given  In  Figs.  3-6.  We  present  records  of  tem¬ 
perature#  pressure#  and  density  for  six  experiments#  and  we  compare  the 
experimental  records  of  density  with  the  results  of  theoretical  calcula¬ 
tions.  Equation  (27)  was  used  for  fractional  densities  up  to  0.90# 
while  equation  (30)  was  used  for  higher  densities.  The  abrupt  change  In 
equations  corresponds  to  an  abrupt  change  In  the  geometry  of  the  model# 
from  highly  connected  open  porosity  to  closed  porosity#  and  a  marked 
discontinuity  in  the  densl flcatlon  rate  might  be  expected.  In  fact#  the 
discontinuity  is  small,  and,  as  shown  In  Figs.  3-6,  It  Is  usually  quite 
difficult  to  see  In  plotted  results.  The  small  size  of  the  discon¬ 
tinuity  suggests  that  there  Is  little  to  be  gained  by  modeling  an  Inter¬ 
mediate  stage  of  denslflcatlon  In  which  "cylindrical  pores"  predominate. 
Neglecting  the  intermediate  stage  is  In  accord  with  the  practice  of 
Swlnkels#  Wilkinson#  Arzt#  and  Ashby  C4]#  who  have  argued  against  the 
use  of  an  Intermediate  stage  model. 

The  Integration  was  carried  out  numerically#  using  an  adaptive 
third-order  Rur.ge-Kutta  method  [16].  From  tests  with  varying  step 
sizes#  It  Is  estimated  that  the  absolute  global  error  In  the  theoretical 
value  of  D  due  to  the  integration  Is  never  larger  than  0.0001#  which  Is 
much  smaller  than  the  expected  error  In  experimental  values  of  D. 

The  only  adjustable  parameters  were  the  initial  value  of  the 
density  and  the  kinetic  factor  f.  For  each  experiment?  a  value  of  f  was 


X-t  z^TuX.t  zrv  r,-v <o-  »  *>.  ar-s; -^w>v.«*^ j 
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chosen  to  provide  a  good  fit  between  theoretical  and  experimental  den¬ 
sity  histories#  while  the  Initial  fractional  density  was  taken  as  the 
value  measured  by  Immersion#  which  is  not  exactly  equal  to  the  density 
given  by  the  dllatometer#  due  to  experimental  error.  In  some  cases#  the 
Immersion  density  was  smaller  than  that  of  a  random  dense  packing  of 
spheres#  Dq,  and  numerical  Integration  starting  at  the  Immersion  density 
was  not  be  possible  due  to  a  singularity  In  equation  (27)  at  D  =  Dq.  In 
those  cases#  an  Initial  fractional  density  of  0.641  was  used. 

The  exact  the  value  of  the  Initial  density  used  In  integration 
actually  has  a  only  very  small  effect  on  the  calculated  density  record# 
since#  In  the  low  density  range#  the  densl flection  rate  drops  rapidly  as 
density  Increases.  As  an  Illustration  of  this#  consider  two  samples 
that  are  being  HIPped  under  Identical#  unchanging  conditions#  and  assume 
that  the  two  samples  are  Identical  except  that  the  Initial  density  of 
one  sample  Is  higher  than  that  c:  t-he  other.  It  will  take  a  certain 
amount  of  time  for  the  less  dense  sample  to  reach  the  Initial  density  of 
the  denser  sample.  Since  the  samples  are  denslfylng  under  Identical 
conditions#  the  densities  of  the  samples  will  be  the  same  function  of 
time  except  for  the  time  offset  just  mentioned.  Although  the  time 
offset  remains  constant#  the  denslflcatlon  rate  drops  so  that  the 
difference  between  the  densities  of  the  two  samples  decreases  as  time 
goes  on.  As  a  qualitative  Illustration  of  the  effect  of  changing  the 
initial  density#  consider  the  results  shown  In  Fig.  7.  The  calculated 
density  history  of  Fig.  3  Is  replotted#  along  with  the  results  of  a  cal¬ 
culation  that  differed  only  In  the  value  of  the  initial  density.  It  Is 
clear  that  the  rather  large  Initial  difference  In  density  quickly 
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becomes  negligible.  Note  also  that  since  the  temperature  and  pressure 
were  not  constant  In  this  calculation,  so  the  time  offset  between  the 
two  curves  Is  not  constant. 


6.  DISCUSSION 


In  the  preceding  discussion,  the  values  of  f  used  In  fitting 
the  theoretical  equations  to  the  data  were  presented  without  any  justif¬ 
ication  that  the  values  are  reasonable.  To  show  that  the  values  are 
Indeed  reasonable,  let  us  estimate  the  values  of  the  variables  that 
comprise  f.  From  scanning  electron  photomicrographs  C17],  It  Is 
estimated  that  R  =  C.35x10*"^m.  Following  Arzt,  Ashby,  and  Verrall  [33, 
we  may  take  bb/b(,  =  V2,  Cj  =  0.5,  P  =  10,  and  bfc  =  1.59*10”i0m,  where 
the  value  of  b^  is  one  third  of  the  Burgers  vector  for  slip  on  the 

{0001}  <1120>  basal  system.  A  reasonable  estimate  of  the  effective 

-79  3 

atomic  volume  of  the  solute  Is  n  =  2.15x10  m  »  which  Is  half  the 

volume  of  the  formula  unit.  Although  the  solute  that  controlled  densif- 

Icatlon  was  not  identified,  chemical  analysis  of  the  powder  revealed 

several  Impurities  for  wh,ch  it  Is  reasonable  to  use  an  atom  fraction  of 

-r  o 

solute  of  C  =  3.x10  The  value  of  D  Is  unknown,  since  It  Is  not 
o  s 

known  what  solute  might  be  responsible  for  controlling  the  motion  of 
grain-boundary  dislocations.  However,  by  combining  the  estimates  above 
with  f  =  3.x104  KPa“^s“*,  a  typlca1  value  obtained  by  curve-fitting,  one 
obtains  a  value  of  D°  =  9.xltf^m^/s.  This  Is  within  the  range  of 
values  of  preexponential  factors  for  solute  diffusion  In  alumina  ?s  com- 
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piled  by  Freer  [18],  so  we  may  conclude  that  the  values  of  f  are  also 
reasonable c 

From  the  captlcns  of  Figs.  3-6,  It  will  be  observed  that  the 
values  of  f  used  In  fitting  theory  to  experiment  vary  by  a  factor  of 
seven.  As  In  Illustration  of  the  effect  of  a  large  change  In  f.  In  Fig. 
8  we  have  replotted  the  theoretical  curve  of  Fig.  5  along  with  a  similar 
curve  for  a  value  of  f  five  times  as  large.  The  two  curves  are  markedly 
different. 

Since  the  variation  In  f  may  seem  to  be  larger  than  might  have 
been  expected.  It  Is  appropriate  to  discuss  possible  causes  of  variation 
In  f.  There  Is  no  evident  correlation  between  the  value  of  f  and  the 
experimental  conditions,  so  we  have  attributed  the  variation  In  f  to 
differences  between  samples.  The  two  factors  In  f  that  would  appear  to 
be  subject  to  the  greatest  variation  are  C  and  D°.  Since  the  powder 

U  P 

was  blended  before  being  formed  Into  samples.  It  Is  not  expected  that 
sampl e-to-sampl e  variation  of  CQ  would  explain  all  of  the  variation  In 
f.  However,  variations  In  D°  are  expected  to  be  much  larger,  since 
solute  dlffuslvltles  In  ceramics  often  depend  very  sensitively  on  the 
concentrations  of  other  Impurities.  Besides  CQ  and  D°»  large  variations 
In  f  may  also  arise  from  relatively  small  variations  In  Eg.  In  our  cal¬ 
culations,  we  used  a  single  activation  energy  for  all  samples.  If  the 
actual  value  of  Eg  for  the  sample  differed  from  the  value  used  In  the 
calculation,  the  Inaccuracy  would  show  up,  greatly  enlarged.  In 
exp(Es/RT).  However,  since  f  was  the  only  kinetic  fitting  parameter, 
changes  in  exp(Es/RT)  would  appear  as  changes  In  f.  The  effects  of  a 
change  In  activation  energy  can  be  quite  significant;  a  variation  of  9% 
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In  Eg  (with  D°  held  constant)  would  not  be  considered  large#  but  It 
would  be  sufficient  to  explain  the  entire  range  of  values  of  f.  There¬ 
fore,  we  conclude  that  sample-to-sample  variations  are  sufficient  to 
explain  the  variation  on  f. 

In  Figs.  3-6,  It  will  be  observed  that  the  theoretical  and 
experimental  curves  are  generally  In  excellent  agreement  up  to  a  frac¬ 
tional  density  of  about  0.9  but  that  the  theory  overestimates  the  den- 
slflcatlon  rate  beyond  this  point.  One  possible  explanation  for  this 
fact  Is  that  the  mechanism  of  denslficatlon  changes  at  this  density. 

But  the  denslficatlon  mechanism  is  determined  primarily  by  the  tempera¬ 
ture  and  tho  driving  force,  and  at  high  densities  the  driving  force  Is  a 
weak  function  of  density.  Therefore,  It  Is  unlikely  that  the  mechanism 
would  change,  under  constant  temperature  and  applied  pressure.  A  more 
plausible  but  less  welcome  explanation  Is  that  the  pore-and-shel 1  model 
described  above  does  not  adequately  describe  the  final  stage  of  donslfl- 
catlon.  This  hypothesis  is  supported  by  the  fact  that  similar  overesti¬ 
mates  of  the  denslficatlon  rate  occur  In  pore-and-shel!  models  of  the 
final  stage  for  other  mechanisms.  Additional  work  on  this  question  Is 
In  progress. 

While  the  agreement  between  theory  and  experiment  Is  generally 
excellent  except  at  high  density,  It  Is  appropriate  to  ask  If  some  other 
denslficatlon  mechanism  could  also  explain  the  data.  Considering  the 
low  temperature  and  small  grain  size  of  this  material,  it  might  be  sup¬ 
posed  that  a  simple  grain-boundary  diffusion  mechanism  might  explain  the 
data.  However,  this  mechanism  has  a  pressure  exponent  of  1:  the  denslfl- 
cation  rate  Is  linearly  dependent  on  the  driving  force  for 
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sure#  such  as  that  presented  In  Figs.  5  and  6#  shows  a  pressure  exponent 
of  about  2  [1.1.  Another  hypothesis  might  be  that  some  lattice  disloca¬ 
tion  mechanism  Is  responsible  for  denslflcatlon#  but  models  of  disloca¬ 
tion  plasticity  generally  give  stress  exponents  of  at  least  3  [19]. 

Also#  most  of  the  experiments  were  performed  at  temperatures  of  1273  to 
1323  K.  It  Is  expected  that  significant  slip  could  occur  only  on  the 
basal  systems  at  these  temperatures  [20]»  and  It  Is  difficult  to  Imagine 
how  significant  denslflcatlon  could  result  with  only  two  active  Indepen¬ 
dent  slip  systems.  Furthermore#  the  activation  energy  for  power-law 
denslflcatlon#  which  Is  generally  thought  to  have  a  dislocation  mechan¬ 
ism#  Is  477  kJ/mol  [21]#  while  the  observed  activation  energy  was  290 
kJ/mol.  This  discrepancy  Is  too  large  to  be  explained  by  experimental 
error.  A  third  possible  explanation  of  the  data  Is  that  two  (or  more) 
simultaneous  mechanisms  are  at  work#  and  that  the  observed  pressure 
exponent  of  2  is  the  result  of  two  mechanisms  with  pressure  exponents  of 
(say)  1  and  3.  Such  an  argument  would  be  difficult  to  defend.  Theoret¬ 
ical  models  of  denslflcatlon  that  give  rise  to  pressure  exponents  of  3 
or  more  generally  Invoke  lattice  dislocation  mechanisms#  and#  as  argued 
above#  lattice  dislocations  can  contribute  little  to  denslflcatlon  at 
the  temperatures  used  In  the  experiments  since  only  basal  slip  can 
occur.  It  is  therefore  concluded  that  a  model  of  denslflcatlon  by 
Interface-reaction  controlled  grain  boundary  diffusion  provides  the  most 
plausible  explanation  of  the  data. 
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Table  of  symbols 


a  average  neck  area 

bfa  Burgers  vector  of  grain  boundary  dislocation 

bn  component  of  Burgers  vector  of  grain  boundary  dislocation 
normal  to  grain  boundary 

c  derivative  of  coordination  number  with  respect  to  R'/R,  c  =  15.5 
Cq  atom  fraction  of  solute  In  bulk 

constant*  =  0.5 
D  fractional  density 

Dq  fractional  density  of  random  dense  packing,  DQ  =  0.64 
D°  dlffuslvlty  of  solute  at  Infinite  temperature 

D  dlffuslvlty  of  solute 

Eg  activation  energy  for  diffusion  of  solute 
F  force  per  unit  length  on  dislocation 
G  shear  modulus 

k  Boltzmann's  constant 

M  dislocation  mobility 

n  constant  for  neck  area  In  Initial  stage,  n  =  11.4 

n^  effective  number  of  grain  boundaries  meeting  at  a  pore 

p  external  pressure 

p#  driving  force  for  denslflcatlon 

Ps  contact  pressure  due  to  surface  tension 

Pth  threshold  contact  pressure  for  dislocation  motion 

R  particle  radius 

R'  current  particle  radius  In  Arzt's  growing-sphere  description 
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Table  of  symbols  (continued) 

Rg  shell  radius 

R  gas  constant 

r  pore  radius 

T  temperature 

t  time 

v  dislocation  velocity 

w  neck  axial  radius  of  curvature 

x  neck  radius 

y  half  particle  center-to-center  distance 

Z  coordination  number 

ZQ  coordination  number  at  DQ,  ZQ  =  7.3 

a  constant,  a  =  0.2 

o  o 

p  ratio  of  solute  concentration  near  dislocation  to  solute 
concentration  In  bulk 
y  surface  tension 

P  dislocation  density 

0  effective  atomic  volume  of  solute 


Fig.  1.  Neck  geometry.  The  two  particles  are  spherical »  while  the  neck 
Is  part  of  a  tangent  circular  torus. 

Fig.  2.  Comparison  of  two  treatments  of  neck  geometry.  Solid  curve  Is 
for  the  toroidal  neck  geometry  of  Fig.  1  with  strict  conservation  of 
mass;  dashed  curve  is  for  the  approximation  given  in  equation  (5).  (a) 

Relative  neck  radius  (x/R)  as  a  function  of  relative  particle  separation 
(y/R).  (b)  Relative  neck  axial  radius  of  curvature  (w/R)  as  a  function 

of  relative  particle  separation  (y/R). 

Fig.  3.  (a)  Pressure  (solid  curve)  and  temperature  (dashed  curve)  his¬ 
tories.  (b)  Comparison  of  experimentally  determined  density  (solid 
curve)  with  Integration  of  kinetic  equations  for  densiflcation  (dashed 
curve)  with  f  =  7.x104  KPa”1s*’1. 

Fig.  4.  (a)  Pressure  (solid  curve)  and  temperature  (dashed  curve)  his¬ 
tories.  (b)  Comparison  of  experimentally  determined  density  (solid 
curve)  with  Integration  of  kinetic  equations  for  densiflcation  (dashed 
curve)  with  f  =  2.8*104  KPa”*s“*. 

Fig.  5.  (a)  Pressure  (solid  curve)  and  temperature  (dashed  curve)  his¬ 
tories.  (b)  Comparison  of  experimentally  determined  density  (solid 
curve)  with  Integration  of  kinetic  equations  for  densiflcation  (dashed 
curve)  with  f  =  l.*104  KPa"1s"1. 

Fig.  6.  (a)  Pressure  (solid  curve)  and  temperature  (dashed  curve)  his¬ 
tories.  (b)  Comparison  of  experimentally  determined  density  (solid 
curve)  with  Integration  of  kinetic  equations  for  densiflcation  (dashed 
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curve)  with  f  =  1.5*104  KPa'V1. 

Fig.  7.  Effect  of  Initial  density  on  denslflcatlon  history.  Temperature 
and  pressure  data  from  Fig.  3  were  Integrated  with  f  =  7.*104  KPa'^s-^, 
but  with  Initial  densities  of  0.641  (solid  curve)  and  0.7  (dashed 
curve). 

Fig.  8.  Effect  of  f  on  denslflcatlon  history.  Temperature  and  pressure 
data  from  Fig.  5  were  Integrated  with  f  =  l.*104  KPa“^s""^  (solid  curve) 
and  f  =  5.*104  KPa'^s-1  (dashed  curve). 
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COMPUTER  SIMULATION  OF  EFFECTS  OF  THE  PORE  SI7C  DISTRIBUTION 
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Battel le»  505  King  Avenue*  Columbus,  Ohio  43201-2693,  USA 

Mstoact 

Most  theoretical  treatments  of  pressure-assisted  densifi- 
cation  of  porous  solids  assume  a  single  size  for  all  pores.  We 
remove  this  this  assumption  and  consider  a  distribution  of  pore 
sizes.  Dissolution  of  intragranular  pores  by  volume  diffusion  and 
dissolution  of  intergranular  pores  by  grain-boundary  diffusion  are 
both  treated.  The  evolution  with  time  of  pore  size  distributions  is 
calculated  for  distributions  that  are  initially  described  by  log¬ 
normal  and  Welbull  functions,  and  differences  in  predicted  behaviors 
are  discussed.  The  pore  size  distribution  Is  then  related  to  two 
Important  quantities:  porosity  and  number  of  pores  per  unit  volume. 
The  assumption  of  a  distribution  of  pore  sizes  Is  found  to  avoid 
certain  unrealistic  predictions  obtained  from  models  with  a  single 
pore  size  such  as  abrupt  disappearance  of  all  pores  and  rapid 
approach  to  full  density. 
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-Introduction 

It  Is  frequently  observed  In  studies  of  denslflcatlon  that 
full  density  Is  difficult  to  achieve.  This  Is  In  direct  conflict 
with  predictions  of  rapid  and  complete  denslflcatlon  derived  from 
theoretical  models  of  denslflcatlon  by  diffusion.  In  this  paper,  we 
show  that  the  models  can  be  brought  into  agreement  with  experience 
by  consideration  of  the  effects  of  a  distribution  of  pr  e  sizes. 

The  importance  of  size-distribution  effects  in  considera¬ 
tions  of  the  evolution  of  a  discrete  second-phase  species  has  long 
been  recognized.  The  porosity  within  a  solid  body  can  be  regarded 
as  a  particular  type  of  second  phase,  and  Its  variation  with  time, 
resulting  from  surface-energy  and  applied-pressure  driving  forces, 
can  be  treated  as  a  type  of  phase  transformation.  In  this  case  as 
well,  size-distribution  effects  can  play  an  Important  role,  for 
example  In  considerations  of  the  manner  in  which  the  overall  volume 
fraction  of  porosity  varies  with  time. 

Presented  below  are  the  results  of  a  modeling  study  of  the 
dissolution  kinetics  of  a  distribution  of  discrete  (l.e.,  non¬ 
overlapping),  spherical  pores  contained  within  a  solid  body.  Both 
intragranular  and  Intergranular  pores  are  considered,  using  well- 
established  expressions  for  the  size-dependent  rate  of  pore  dissolu¬ 
tion  resulting  from  volume  diffusion  and  grain-boundary  diffusion, 
respectively,  of  vacancies  away  from  the  pore  surface.  The  evolu¬ 
tion  with  time  of  the  pore  size  distribution  Is  calculated,  starting 
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from  various  assumed  initial  distributions*  and  from  this  the 
corresponding  variation  of  the  overall  porosity  of  the  solid  is  cal¬ 
culated  and  related  to  size-distribution  effects. 

Because  of  the  fact  that  both  pore-surface-energy  and 
applied-pressure  driving  forces  are  considered*  the  problem  Is  not 
analytically  tractable  for  either  Intragranular  or  Intergranular 
pores.  Consequently*  solutions  are  generated  numerically. 

It  is  clear  that  the  assumption  of  a  system  of  discrete 
pores  limits  our  consideration  to  a  solid  for  which  the  amount  of 
porosity  Is  relatively  low.  This  corresponds  to  what  Is  commonly 
referred  to  as  the  "final  stage"  of  denslflcatlon.  A  quantitative 
evaluation  of  the  probability  that  a  given  pore  is  actually  discrete 
can  be  carried  out,  as,  for  example*  has  been  done  Cl]  for  the  case 
of  a  uniform  size  distribution  of  spheres  distributed  at  random 
within  three-dimensional  space. 

J^or e-D Issolut  1  on  Models 

Over  the  years*  many  investigators  have  developed  models 
for  the  growth  and  dissolution  of  pores  contained  within  a  solid 
body.  (Actually,  models  for  pore  growth  can  be  used  to  describe 
dissolution  as  well*  If  appropriate  modification  Is  made  of  the 
applied-stress  term  In  the  growth-rate  expression,  as  has  been  dis¬ 
cussed  by  Greenwood  [23. )  In  addition,  a  variety  of  kinetic  mechan¬ 
isms  has  been  considered.  For  purposes  of  the  present  analysis* 
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relatively  simple  models  for  pore  dissolution  are  used  which 
nevertheless  serve  well  to  Illustrate  the  manner  In  which  size- 
distribution  effects  can  influence  the  overall  kinetics  of  porosity 
reduction. 


-Intr  <&r  nJEaEfis. 

For  the  volume-dl f fusion-controlled  dissolution  of  an  Iso¬ 
lated,  spherical*  Intragranular  pore  (1.e.»  neighboring  pores  are 
assumed  not  to  Influence  one  another)*  one  can  derive  an  expression 
for  the  rate  of  pore  dissolution  by  assuming  that  the  vacancy- 
concentration  field  within  the  solid  outside  the  pore  satisfies  the 
Laplace  equation.  One  thus  obtains  (e.g.,  [3,4]), 


or> 

dR=_— ^fp  +  211  {1) 

dt  kTR  R  J 

where  R  Is  the  pore  radius  at  time  t»  Q*  Dy,  and  r  are  the  atomic 
volume,  the  volume  sel f-dlffusivlty,  and  the  surface  energy  of  the 
solid,  P  Is  the  externally  applied  pressure,  k  is  Boltzmann’s  con¬ 
stant,  and  T  Is  the  absolute  temperature.  Three  assumptions 
Inherent  in  Equation  1  are  (a)  that  the  vacancy  supersaturation  In 
the  lattice  Is  zero,  lb)  that  the  "effective"  applied  pressure 
inside  the  solid  is  unaffected  by  the  presence  of  the  porosity,  and 
(c)  that  no  gases  exist  Inside  the  pores.  Relaxation  of  assumption 
(a)  has  been  considered  by  Geguzln  and  Llfshitz  [3]  and  others; 
relaxation  of  (b)  and  (c)  has  been  discussed  by  Markworth  [4]  and 


others 
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In  order  to  simplify  the  analysis  which  follows,  we  re' 
express  Equation  1  In  terms  of  dimensionless  variables,  l.e.,  a 
pore-size  parameter  p  and  a  time  parameter  %  defined  as 


In  terms  of  which  Equation  1  assumes  the  much  simpler  form 


=  _i  r 

dt  p  L 


1  + 


(2) 

(3) 

(4) 


thickness  of  the  grain  boundary*  and  all  other  symbols  are  as 
defined  for  Equation  1.  Again,  It  Is  assumed  that  the  ••effective" 
applied  pressure  inside  the  solid  is  unaffected  by  the  presence  of 
porosity  and  that  no  gases  exist  inside  the  pores.  Effects  of  gases 
Inside  the  pores  could  be  considered  [5]  and  are  analogous  the 
corresponding  effects  for  Intragranular  pores. 

It  Is  again  convenient  to  express  Equation  5  in  terms  of 
dimensionless  variables.  We  use  the  same  size  and  time  parameters, 
p  and  t,  given  by  Equations  2  and  3,  but  we  define  a  new  Himens1on- 
less  factor,  a,  as 

a  =  (6) 

and  obtain  the  following  simpler  form  for  Equation  5: 


The  rate  of  pore  shrinkage  -dp/dt  as  a  function  of  pore  radius  p  Is 
plotted  In  Fig.  1  for  both  Intergranular  and  Intragranular  pores,  as 
calculated  from  Equations  4  and  7  with  a  value  of  a  =  1.  From  Fig. 
1,  It  can  be  seen  that  the  rate  of  pore  shrinkage  decreases  as  pore 
radius  increases,  and  that  the  dependence  of  shrinkage  rate  on  pore 
size  Is  stronger  for  Intergranular  pores  than  it  is  for  Intragranu¬ 
lar  pores.  Note  that  the  shrinkage  rate  for  Intergranular  pores  Is 
dependent  on  a,  and  that  the  corresponding  curve  In  Fig.  1  would  be 
raised  or  lowered  If  a  different  value  of  a  hc.d  been  chosen. 


Clearly*  pore-dissolution  models  other  than  those 
represented  by  Equations  4  and  7  could  have  been  used.  These  two 
particular  models  were  chosen  because  they  are  particularly  amenable 
to  the  size-distribution  analyses  developed  below  while  still  con¬ 
taining  a  description  tf  some  essential  features  of  the  physics  of 
the  respective  dissolution  processes. 

One  Interesting  case  not  mentioned  above  Is  the  oft-quoted 
model  of  Hull  and  Rlmmer  [6]  for  the  stress-induced  growth  of 
grain-boundary  voids.  Using  their  model,  and  Including  only 
applied-pressure  and  surface-energy  effects  as  arlvlng  forces  for 
pore  dissolution,  one  obtains  an  expression  equivalent  to  Equation  4 
as  the  dimensionless  representation,  noting  that  Equation  4  was 
derived  for  a  model  of  Intragranular  pores.  The  only  difference 
would  be  the  replacement  of  Dy  in  Equation  1  with  Db6/(2a)  where  a 
Is  the  mean  separation  between  pores.  Consequently,  the  analysis  of 
Intragranular  porosity  based  on  Equation  1  is  mathematically 
equivalent  to  that  of  Intergran'^^r  porosity  based  on  the  Hull- 
Rlmmer  model  [6].  It  should  d,  however,  that  Important 

corrections  to  the  Hull-Rimmer  inr-  have  been  made,  e.g.,  by  Weert- 
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3.  Size-Distribution  .Kinetics 

3.1.  Evolution  of -the  Pore  Size  Distribution 

We  define  the  size-distribution  function  for  a  system  of 
pores  as  f(R,  t),  such  that  f(R,  t)dR  Is  the  concentration  of  pores* 
at  time  t>  having  radius  within  the  range  R  to  R  +  dR.  Clearly,  It 
Is  assumed  In  this  definition  that  the  size  distribution  exhibits  no 
spatial  dependence,  on  the  average.  Assuming  that  no  ••sources"  or 
"sinks"  for  pores  exist,  one  can  show  that  the  function  f  satisfies 
the  following  continuity  relation: 

ft  +  3R<vfl  *  0  (8) 

where  v  is  a  function  that  describes  the  rate  at  which  a  given  pore, 

which  can  be  considered  as  a  point  existing  in  a  one-dimensional 
pore-radius  space,  moves  through  that  space.  (Obviously,  the  defin¬ 
ition  of  v  is  equivalent  to  dR/dt  which  Is  given  In  Equations  1  and 
5  for  the  two  models  under  consideration  here.)  Expressing  Equation 
8  in  terms  of  the  dimensionless  parameters  p  and  x, 

J*  +  ~(vf)  =  0  (9) 

OX  f*p 

where  now  f(p,  t)dp  Is  the  concentration  of  pores,  at  time  x  having 
rad’us  within  the  range  p  to  p  +  dp,  and  v  Is  now  equivalent  to  the 
"velocity",  dp/dt.  In  p-space. 

The  "flux"  of  pores  crossing  a  given  point  of  the  one¬ 
dimensional  p-space  Is,  In  general,  equal  to  v(p»  r)f(p,  t).  For 
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the  pore-dissolution  models  under  consideration  here»  v  has  explicit 
dependence  only  upon  p  (see  Equations  4  and  7). 

Let  us  examine  the  behavior  of  the  pore  flux  with  respect 
to  a  moving  coordinate  frame  In  our  one-dimensional  p-space. 
Specifically*  let  p'  be  a  coordinate  which  moves  along  the  p-axls  at 
the  same  rate  as  a  pore  having  Instantaneous  size  p.  Hence*  p’  Is 
given  through  the  relation 


(10) 


where  pQ  is  the  size  of  the  pore  at  some  given  Initial  time  t0  and 
where  we  again  take  v  to  depend  only  upon  p.  From  Equation  10  we 


obtain 


Ilf'  '"p'1 

It  now  follows  that* 


(11) 


[v(p»)f(p»,  t)]  (12) 

♦  vipnrjiL-np.,  t,  i^tJ4,(p.,  *>]  . 

Combining  Equations  11  and  12* 

[v(p»)f(p»,  t)]  (13) 

'  »<*'>  L^dp^  f<c’'  +  «►’’  *>  +-£  ««'■  *>]  • 

Now*  the  continuity  equation.  Equation  9,  must  a’lso  be  satisfied 
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with  respect  to  the  moving  coordinate  system,  in  which  case  the  sum 
of  the  terms  within  brackets  on  the  right  hand  side  of  Equation  13 
Is  zero.  Consequently, 

[v(p’)f(p»,  t)3  =  0  .  (14) 
An  alternative  way  of  expressing  Equation  14  is  the  following: 


v(p)f(p,  t)  =  v(p)f(pQ,  t)  (15) 

where  p  (dropping  the  prime  symbol)  Is  understood  to  be  the  size  of 
a  pore,  at  time  x  that  had  size  pQ  at  an  earlier  time  f0.  Taking 
T0  =  0  and  g(PQ)  =  f(pQ,  0),  we  obtain 


f(p,  x) 


(16) 


Equation  16  can  be  used  to  relate  the  size  distribution  at  time 
x  =  0  to  that  at  some  late  time  x  >  0.  An  alternative  derivation  of 
Equation  15  is  presented  In  the  appendix. 

In  order  to  apply  Equation  16  to  cases  of  Interest  here, 
one  can  use  Equation  10  (again,  and  from  now  on,  dropping  the  prime 
symbol)  together  with  the  given  expressions  for  v(p).  This  can  be 
done  in  closed  form  for  the  models  under  consideration  here.  In 
particular,  substituting  the  right-hand  side  of  Equation  4  for  v(p) 
into  Equation  10  and  Integrating,  one  obtains  the  following  expres¬ 
sion  for  intragranular  pores: 


o 


P)(P0  +  P  -  2)  +  In 


LUm] 

1  +  p  J  • 


•c 


(17) 


Likewise*  substituting  the  right-hand  side  of  Equation  7  for  v(p) 
Into  Equation  10  and  Integrating,  one  obtains  the  following  expres¬ 
sion  for  Intergranular  pores: 

ax  =  3(po  '  p3)  +  (po  "  p)  “  (po  “  p2>  +  ln  j  *  (18) 

Unfortunately,  for  given  values  of  r  and  pQ,  the  value  of  p  can  be 
calculated  from  either  Equation  17  or  18  only  through  some  numerical 
procedure.  Given  this  fact,  however,  the  problem  of  calculating  the 
evolution  of  the  pore  size  distribution  Is,  In  principle,  solved. 

The  above  analysis  can  be  cast  In  somewhat  simpler  form  by 

taking 

Po  =  P  +  A  (19) 

where  A  >  0  since  the  pores  are  dissolving.  In  this  form.  Equations 
16,  17  and  18,  respectively,  become 

f(p,  t)  =  V-^(-p-)A^ ( P  +  A)  (20) 

x  =  A(P  -  1  +  £a)  +  In  -+pA]  (21) 

ax  =  A(p2  +  pA  +  *A2  +  1  -  p  -  ^A)  -  In  [  1  -f  -+  p  A  ]  •  (22) 

For  given  p  and  x ,  Equations  21  and  22  can  be  solved  numerically  for 
the  corresponding  values  of  A  for  Intragranular  and  intergranular 
pores,  respectively. 
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3-c2^  -Ev-Ol  jitlflflL  _Qf  -D  1st  r  1-b  u  t  ton-_Fji.net  ton.  ■Moments 

Properties  of  physical  Interest  are  related  to  moments  of 
the  size-distribution  function,  rather  than  to  the  distribution 
function  per  se>  defining  the  1-th  algebraic  moment  (i.e.,  the  1-th 
moment  about  p  =  0)  as 

03 

M.(x)  =  /p^ffp,  *c ) dp  .  (23) 

1  0 

For  example,  Mg(x)  is  the  net  concentration  of  pores  within  the 
solid,  and  the  volume  fraction  occupied  by  porosity  Is  linearly  pro¬ 
portional  to  M3(t). 

There  are  different  ways  In  which  M^(r)  can  be  calculated 
(see,  e.g.,  [8,  9]).  As  d1ccussed  below,  the  approach  used  here  Is 
to  calculate  f(p,  x)  using  a  general  approach  described  above,  and 
then  numerically  Integrating,  using  Equation  23  as  a  basis,  to 
determine  the  moments  of  interest. 

One  case  is  of  particular  Interest,  namely,  that  Involving 
Mq.  It  can  easily  be  shown  that  an  alternative  expression  for  Mg  Is 

CD 

Mg(x)  =  /  g(p)dp  (24) 

A 

o 

where  the  quantity  Aq  in  Equation  24  Is  the  value  of  A,  obtained 
from  Equation  21  or  22,  corresponding  to  p  =  0.  Thus,  for 
Intragranular  and  Intergranular  pores,  respectively. 


rr^jzjxr^ jr^n.  -"wv raurcura./  rv*. jw»  ruir* Ajwi.'w^mrun  7\p»  fwrJiAiirkA ajuw  .:\h  rd*  tw  r^vTjR  f\j* r^/^L. 
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T  =  A0(^A0  -  1)  +  In (1  +  A0> 


"  ■  4o(34o  -  2Ao  -  11  -  1"(1  +  V  • 


(25) 

(26) 


The  values  of  Mg  can  be  calculated#  for  a  given  value  of  *c»  simply 
by  numerically  determining  the  pertinent  value  of  Aq  from  Equation 
25  or  26#  and  then  performing  the  integration  Indicated  in  Equation 
24.  The  result  obtained  by  this  approach  should  yield  the  same 
result  as  that  obtained  from  the  approach  previously  described  based 
on  the  use  of  Equation  23.  In  fact#  comparison  of  the  values  of  Mg 
calculated  using  the  two  different  approaches  would  serve  as  a  use¬ 
ful  check  of  the  accuracy  of  the  numerical  procedures  and  was  Indeed 
used  for  this  purpose  in  the  examples  described  below. 


-4..  jLa_5p_ectflc .  Examp  1  es 

The  methods  of  calculating  size-distribution  kinetics 
described  above  may  be  applied  to  a  wide  variety  of  Initial  pore- 
size  distributions.  We  have  applied  them  to  a  Welbull  distribution 


m 


|M  I 

g(p)  =  mpm“1e"p  .  (25) 

with  m  =  2#  2.5#  and  3#  and  to  a  log-normal  distribution 


g(p>  = 


exp 


rXln.U/h-11' 


(2c2) 


(26) 


V2ncp 

where  we  have  taken  b  =  0.8  and  c  =  0.5.  As  will  be  seen  from  the 


figures  below#  these  distributions  are  generally  similar  in  appear¬ 
ance.  The  most  Important  difference  among  them  is  that  the  log¬ 
normal  distribution  has  a  long  large-radius  tall,  while  the  large- 
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radius  tails  of  the  Welbull  distributions  are  smaller  and  decrease 
with  Increasing  m. 

All  of  these  functions  are  normalized#  that  Is#  they  have 
Mq  =  1.  However#  due  to  the  differences  between  the  functions#  It 
Is  generally  not  possible  to  normalize  the  third  moments  simultane¬ 
ously.  For  the  Welbull  distributions#  we  have  =  T(l  +  3/m)  while 
for  the  log-normal  distribution,  =  fc^exp(9c^/2) . 

In  Figs.  2-4  we  present  plots  of  the  pore  size  distribu¬ 
tion  function  for  selected  times  and  various  Initial  distributions 
as  calculated  using  Equation  4  for  the  kinetics  of  dissolution  of 
Intragranular  pores.  Results  for  the  Welbull  distribution  with 
m  =  2.5  are  not  plotted  but  were  Intermediate  between  the  results 
for  m  =  2  and  m  =  3.  From  Fig.  2#  It  will  be  noted  that  the  posi¬ 
tion  of  the  peak  of  the  distribution  function  moves  toward  larger 
radii  as  time  progresses.  A  similar  but  smaller  effect  can  be  seen 
In  Fig.  3#  while  In  Fig.  4#  position  of  the  peak  clearly  shifts 
toward  smaller  radii.  This  behavior  Is  In  marked  contrast  to  that 
of  standard  models  In  which  all  pores  are  the  same  size  and  all 
pores  shrink  together:  for  denslflcatlon  to  occur#  the  pore  radius 
must  shrink.  If  there  Is  a  distribution  of  pore  slz^s#  however#  the 
peak  of  the  distribution  may  move  In  either  direction#  depending  on 
the  shape  of  the  distribution  and  the  velocity  function.  The  Influ¬ 
ence  of  the  velocity  function  may  be  seen  In  Figs.  5-7,  which  differ 
from  Figs.  2-4  only  In  that  the  velocity  for  intergranular  pores  was 
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used.  For  all  calculations  with  intergranular  pores,  we  took  o  =  1. 
As  seen  In  Figs.  5-7,  the  peak  shifts  strongly  to  higher  radii  for 
the  log-normal  distribution,  somewhat  less  strongly  for  the  Welbull 
distribution  with  m  =  2,  and  only  slightly  for  the  Welbull  distribu¬ 
tion  with  m  =  3.  It  Is  clear  from  the  figures  that,  qualitatively, 
if  the  magnitude  of  the  velocity  function  decreases  strongly  with 
Increasing  radius  and  the  initial  distribution  has  a  long  large- 
radius  tall,  the  peak  of  the  distribution  function  will  shift  toward 
the  right,  since  the  small  pores  disappear  quickly,  leaving  the 
larger  pores  almost  unchanged.  Conversely,  If  the  velocity  function 
depends  weakly  on  radius  and  the  Initial  distribution  Is  narrow,  the 
pores  will  shrink  together  and  the  peak  of  the  distribution  will 
move  to  the  left. 

While  the  pore  size  distribution  functions  would  be  diffi¬ 
cult  to  determine  experimentally,  two  quantities  that  are  more 
easily  measured  are  the  number  of  pores  per  unit  volume  and  the 
porosity.  These  two  quantities  are  proportional  to  and  Mg, 
respectively.  The  moments  are  plotted  as  functions  of  time  for  all 
four  distributions  and  both  velocity  functions  In  Figs.  8-11.  For 
comparison,  we  have  also  plotted  the  third  moments  for  a  distribu¬ 
tion  In  which  all  pores  are  the  same  size.  An  Initial  pore  radius 
of  p  =  < 3/4jt )  1/3  was  usecj# 

If  It  Is  assumed  that  all  pores  have  the  same  size,  all 
pores  vanish  simultaneously,  and  changes  discontlnuously.  A  more 
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realistic  picture  Is  seen  In  Figs.  8  and  10:  the  number  of  pores 
gradually  decreases  toward  zero.  At  times  up  to  r  =  0.2,  the  number 
of  pores  decreases  at  comparable  rates  for  all  four  distributions* 
reflecting  the  general  similarity  of  the  central  portion  of  the  dis¬ 
tributions.  By  t  =  1,  however*  the  distributions  are  dominated  by 
what  was  originally  the  large-radius  tall  of  the  distribution,  and, 
for  both  velocity  functions,  the  log-normal  distribution  has  the 
largest  number  of  pores,  followed  by  the  Welbull  distributions  In 
order  of  increasing  m. 

From  Figs.  9  and  11,  we  note  again  the  unrealistic  results 
obtained  by  assuming  a  single  pore  size:  denslflcatlon  proceeds 
rapidly  to  completion.  At  short  times,  the  rates  of  change  of  the 
third  moment  as  obtained  by  assuming  .  .rlbutlon  of  pore  sizes 
agree  relatively  well  with  those  for  a  single  pore  size  and  with 
each  other.  However,  the  results  soon  begin  to  diverge.  For  a  sin¬ 
gle  size  of  Intragranular  pores,  It  Is  predicted  that  all  porosity 
will  vanish  at  x  ~  0.247,  but  all  of  the  calculations  with  a  distri¬ 
bution  of  sizes  give  significant  remaining  porosity  at  this  time. 

The  contrast  is  even  stronger  In  the  case  of  intergranular  pores. 

The  calculation  with  a  single  pore  size  shows  all  porosity  vanishing 
at  x  -  0.196,  but,  for  the  log-normal  dlrtrlbutlon  and  Welbull  dis¬ 
tribution  with  m  =  2,  more  than  half  of  the  original  porosity 
remains.  For  both  velocity  functions,  the  calculations  on  a  log¬ 
normal  distribution  give  significant  remaining  porosity  at  x  -  1. 
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S.  Discussion 

Size-distribution  effects  associated  with  pores  In  solids 
have  been  studied  by  other  Investigators  as  well.  For  example* 
Tomandl  [10]  considered  pore  shrinkage  resulting  from  sintering 
without  any  externally  applied  pressure.  The  expression  he  used  to 
describe  the  rate  of  pore  dissolution  by  diffusion  along  grain  boun¬ 
daries  was 

#  *  (27) 

d‘  ^ 

where  c0  is  a  parameter  that  Is  dependent  upon  temperature. 

Clearly,  Equation  27  Is  mathematically  equivalent  to  Equation  1  If  P 
is  set  equal  to  zero  In  the  latter.  They  then  differ  only  in  the 
form  of  the  coefficient  of  the  term.  For  this  limiting  case* 
the  time-dependent  size-distribution  function  can  be  derived  in 
closed  form,  as  Tomandl  has  shown  CIO].  In  other  studies  Cll*  12], 
the  evolution  of  a  cavity  size  distribution  under  an  applied  tensile 
stress  was  studied.  The  rate  of  change  of  pore  size  was  here 
described  by  an  expression  equivalent  to  our  Equation  5*  with  P 
replaced  with  -o,  where  a  Is  the  tensile  stress  at  and  perpendicular 
to  the  grain  boundary.  It  was  assumed  that  a  Is  equal  to  the  exter¬ 
nally  applied  tensile  stress,  an  assumption  which  (as  was  pointed 
out  C12] )  is  not  generally  valid.  The  evolution  of  the  cavity  size 
distribution  was  evaluated  numerically  using  an  approach  that  was 
analogous  to  that  used  here. 
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It  Is  a  common  experience  in  sintering  and  hot  isostatic 
pressing  to  observe  that  it  is  difficult  to  achieve  full  density. 
This  Is  In  marked  contrast  to  the  predictions  of  rapid  and  complete 
densiflcatlon  obtained  from  standard  models  derived  under  the 
assumptions  of  dlffuslonal  control  and  a  single  pore  size.  We  have 
shown  that  more  realistic  results  may  be  obtained  by  assuming  that  a 
distribution  of  pore  sizes  exists  In  the  material.  We  have  also 
shown  that  the  most  likely  pore  size  can  increase  even  as  denslfica- 
tion  proceeds  and  all  pores  shrink. 

It  Is  beyond  the  scope  of  this  paper  to  relate  pore  size 
distributions  to  particle  size  distributions  and  particle  packings. 
However,  It  Is  clear  from  this  work  that  behavior  in  the  limit  of 
long  time  (and  high  density)  is  controlled  by  large  pores.  It  Is 
expected  that  these  large  pores  would  be  found  near  large  Initial 
particles  or  fully  denslfled  agglomerates. 
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AiuteMtxiALtern_a.ttv.e- _P_e.rlY.at ton  _of_  Equation  15 

In  this  appendix,  we  present  an  alternative  derivation  of 
Equation  15,  one  that  Is  perhaps  simpler  than  that  developed  in  the 
text,  but  which  nevertheless  is  mathematically  rigorous. 


iru  wi.wjsa  w^nLw^v^  xv  irzxrimurxw^wuttJ^TA?:  vwlwuvw  tmsuwoiTO 


Consider  a  pore  that  has  size  Pq  as  time  tq  and  size  p  >  0 
at  some  later  ilme  t.  The  relationship  between  p  and  x  is  given  by 
Equation  10  (dropping  the  prime  symbol  In  Equation  10).  Likewise, 


consider  another  pore  that  has  size  pQ  +  dPQ  at  time  xq  and  size 


p  4  dp  at  time  t.  For  this  case.  Equation  10  becomes 


1 

P  +  dp  . 

o  “  '  v(  p)  • 

Co  +  dpo 

(A.l) 

sT 

if 

Subtracting  Equation  10  from  A.l, 

& 

U 

en 

0=  p  )  dP  -4<L_  .  /JL 
u  J  v(p)  J  V(P)  * 

po  +  dpo  po 

(A. 2) 

te 

or 

s 

B 

P  +  dp  .  p  . 

f  _d£_  _  f  dp 

o  ho  v(p)  D  v(P)  * 

P0  'r  dPn  Pn 

(A.3) 

v  u  u 

For  Infinitesimally  small  dpQ  and  dP, 

Equation  A.3  can  be  expressed 

R 
}*  ^ 

as 

-djP-.  _ 


v(p)  v(pQ)  * 


(A. 4) 


The  concentration  of  pores  at  time  rQ  within  the  interval 


p0  to  po  +  dpo  ^po'  To)dpo  and  at  t*fne  T  w1th*n  the  Inter¬ 


val  p  to  p  +  dp  is  f(p,  x)o‘p.  Clearly,  these  must  be  equal,  1.e„, 


f(p,  x)dp  =  f(p0,  -co)cJP0  . 


(A. 5) 


Eliminating  dp  and  dp  from  Equation  A. 5  by  application  of  Equation 


A. 4,  we  find  that 


3m  sss  tpsn  m3  @  @es> 
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v(p)f(p,  r)  =  v(po)f(pQ,  t0)  (A. 6) 

and  we  see  that  Equations  A. 6  and  15  are  identical. 
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Fig.  1.  Normalized  shrinkage  rate  as  a  function  of  normalized  pore 
radius  for  two  types  of  pores.  Solid  line:  tntragranular  pores. 
Dashed  line:  Intergranular  pores. 

Fig.  2.  Pore  size  distribution  function  for  int. agranular  pores  at 
several  times.  Initial  distribution  Is  log-normal. 

Fig.  3.  Pore  size  distribution  function  for  Intragranular  pores  at 
several  times.  Initial  distribution  Is  Welbull  with  m  =  2. 

Fig.  4.  Pore  size  distribution  function  for  intragranular  pores  at 
several  times.  Initial  distribution  Is  Welbull  with  m  =  3. 

Fig.  5.  Pore  size  distribution  function  for  intergranular  pores  at 
several  times.  Initial  distribution  Is  log-normal. 

Fig.  6.  Pore  size  distribution  function  for  intergranular  pores  at 
several  times.  Initial  distribution  Is  Welbull  with  m  =  2. 

Fig.  7.  Pore  size  distribution  function  for  intergranular  pores  at 
several  times.  Initial  distribution  Is  Welbull  with  m  =  3. 

Fig.  8.  Zeroth  moment  of  the  pore  distribution  function  (propor¬ 
tional  to  number  of  pores)  as  a  function  of  time  for  intragranular 
pores  with  four  Initial  pore  size  distributions.  Solid  line:  log¬ 
normal.  Short-dashed  line:  Welbull ,  m  =  2.  Dot-dashed  line: 
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Weibull,  m  =  2.5.  Long-dashed  lines  Weibull,  m  =  3. 

Fig.  9.  Third  moment  of  the  pore  distribution  function  (proportional 
to  porosity)  as  a  function  of  time  for  Intragranular  pores  with  five 
Initial  pore  size  distributions.  Solid  line:  log-normal.  Short- 
dashed  line:  Weibull,  m  =  2.  Dot-dashed  line:  Weibull,  m  =  2.5. 
Long-dashed  line:  Weibull*  m  =  3.  Dotted  line:  single  pore  size. 

Fig.  10.  Zeroth  moment  of  the  pore  distribution  function  (propor¬ 
tional  to  number  of  pores)  as  a  function  of  time  for  Intragranular 
pores  with  four  initial  pore  size  distributions.  Solid  line:  log¬ 
normal.  Short-dashed  line:  Weibull,  m  =  2.  Dot-dashed  line: 
Weibull,  m  =  2.5.  Long-dashed  line:  Weibull,  m  =  3. 

Fig.  11.  Third  moment  of  the  pore  distribution  function  (propor¬ 
tional  to  porosity)  as  a  function  of  time  for  Intragranular  pores 
with  five  Initial  pore  size  distributions.  Solid  line:  log-normal. 
Short-dashed  line:  Weibull,  m  =  2.  Dot-dashed  line:  Weibull,  m  = 
2.5.  Long-dashed  line:  Weibull,  m  =  3.  Dotted  line:  single  pore 
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